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I . Introduction 


In  July  1977  the  Naval  Postgraduate  School  (NFS)  conducted  a cruise 
aboeurd  the  R/V  Acania  of  approximately  two  weeks  duration  (7/17  to  7/28) 
to  the  lios  Tbigeles  air  basin.  The  basin  is  of  great  Interest  for  meteor- 
ological studies  because  of  its  configuration,  which,  in  conjunction  with 
the  prevalent  marine  inversion  causes  a severe  air  pollution  problem. 

The  purpose  of  this  cruise  was  to  conduct  initial  overwater  studies  of 
the  area  within  40  nautical  miles  of  shore  from  Santa  Barbeura  to  Long 
Beach.  There  is  a need  for  this  type  of  work  since  very  little  data  has 
been  collected  on  the  seaward  boundary  of  the  air  basin.  This  mecuis  that 
there  is  insufficient  data  available  to  establish  boundary  conditions  for 
current  air  pollution  models.  Since  the  westward  boundary  of  the  area  is 
feir  out  to  sea,  data  must  be  collected  from  some  type  of  overwater  plat- 
form. 

One  of  the  main  pxirposes  of  this  cruise  was  to  collect  turbulence  data 
on  the  boundaries  of  the  California  Air  Resources  Board  (CARB)  numerical 
air  pollution  diffusion  model  zmd  at  several  locations  within  the  boundary. 
This  data  was  to  be  collected  at  several  different  times  of  day,  and  if 
possible,  under  different  meteorological  conditions.  At  the  same  time  a 
fairly  con^)lete  set  of  parameters  were  measured  in  order  to  specify  the 
local  atmospheric  conditions.  These  data  will  be  used  to  adjust  the 
peurameterization  of  the  CARB  model. 

The  second  main  endeavor  of  the  NPS  effort  was  to  participate  in  an  air 
flow  trajectory  e]q>eriment  conducted  by  the  California  Institute  of  Tech- 
nology (CIT) . In  this  experiment  a tracer  gas  (SFg)  was  released  from  the 
smoke  stacks  of  the  El  Segundo  power  plant  emd  the  air  flow  was  traced  by 
monitoring  stations  along  the  shore  and  by  the  R/V  Acania.  Thus,  the  ship 


provided  a portion  of  the  tracer  effort  and  the  overwater  meteorological 
support.  Of  course,  the  meteorological  data  was  also  used  to  expcuid  the 
data  base  for  the  model  studies. 

A complete  list  of  the  various  activities  performed  on  the  cruise  is 
given  in  Section  II. 

This  report  will  describe  in  detail,  and  report  the  results  for  the 
NPS  effort  on  this  cruise.  Other  groups,  in  particulsu:  the  Air  Resources 
Board,  were  major  participiuits  in  the  cruise  and  their  results  are  reported 
elsewhere.  Here  we  describe  the  other  projects  only  very  briefly,  and  no 
results  are  reported.  Of  course,  the  R/V  Acania  was  the  vehicle  for  all 
overwater  experiments. 

This  was  planned  as  the  first  of  a number  of  cruises  whose  purpose 
would  be  to  collect  overwater  data  which  would  directly  relate  to  air  pol- 
lution studies.  The  period  of  time  chosen  was  not  necessarily  the  best 
for  performing  some  of  the  studies  (in  paurticular  the  release  of  SFg  from 
the  power  plant)  but  was  a target  of  opportunity.  However,  since  we  hope 
to  conduct  several  such  cruises  during  different  times  of  the  year  in 
order  to  test  a wide  range  of  conditions,  the  dates  used  are  entirely 
appropriate.  Also,  as  shall  be  seen  below,  the  prevailing  conditions  were 
excellent  for  all  of  the  e^erlments  planned. 
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II.  Participants,  Ej^eriments  Performed 


The  various  organizations  other  than  NPS  which  participated  in  this 
cruise  were: 

California  Air  Resources  Boeurd  (GARB) 

Calspan  Corporation 

California  Institute  of  Technology  (CIT) 

Naval  Weather  Service  Facility  (NWS) 

Science  Application,  Inc.  (SAI) 

Western  Oil  & Gas  Association  (WOGA) 

The  principal  activities  of  these  groups  were  as  follows: 

CARS:  Collect  data  on  air  pollutants,  coordination  between 

shipboard  and  shore  personnel,  logistic  support  in 
the  Southern  California  area. 

Calspan:  Collect  aerosol  data,  maintain  meteorological  conditions 

log. 

CIT:  Conduct  release  ^md  collection  of  tracer  gases  both  on 

ship  zmd  shore. 

NWS:  Perform  radiosonde  releases  twice  daily. 

SAI:  Observe  shipboard  operations  2uid  maintain  a log  for  WOGA. 

WOGA:  Coordinate  with  oil  tankers  and  drilling  platforms. 

The  various  experiments  and  the  dates  and  times  on  which  they  were 
performed  were: 

17/0900  - 18/0230  Monitor  open  ocean  aerosol  and  pollutemt  back- 
ground levels. 

18/0300  - 18/1000  Monitor  aerosol  and  pollutwts  in  Semta  Barbara 

Cheumel,  measure  oil  platform  emission  levels. 
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18/1945  - 19/0345 
T1  19/1100  - 19/1200 


19/1430  - 19/1500 

Ml 

19/1545  - 19/1730 

T2 

20/0430  - 20/0820 

P 

20/1710  - 20/1800 

M2 

20/1820  - 20/2200 

M3 

21/0000  - 21/2115 

G1 

21/2300  - 22/1100 

T3 

22/1700  - 22/2100 

22/2100  - 23/0830 

S 

23/1400  - 23/1800 

G2 

23/2300  - 24/1030 

25/2100  - 26/0320 

G3 

26/0530  - 26/1730 

27/0000  - 27/0445 

27/2030  - 27/2230 

28/0500  - 28/0715 

Near  shore  turbulence  sttidies  In  Santa  Monica  Bay. 
Monitor  tanker  emissions  during  non-transfer 
period. 

Cross  check  ship  amu  shore  instruments  at  Catalina. 
At  sea  turbulence  data. 

Monitor  emissions  during  tanker  treuisfer  operation. 
Monitor  drilling  platforms. 

Turbulence  data  on  Southern  edge  of  CABB  model. 
Santa  Monica  Bay  CABB  model  studies. 

Track  tracer  gas  released  from  power  plant  on 
shore . 

Monitor  emissions  during  tanker  transfer  operation. 
Monitor  pollutants  and  aerosols  along  coast  from 
Del  Mar  area  to  Long  Beach  2u:ea. 

Check  turbulence  results  for  various  ship  maneuvers 
Track  tracer  gas  released  from  power  plant  on  shore 
Turbulence  study  during  light  wind  conditions. 
Belease  of  tracer  gas  from  ship  along  shipping 
lanes  west  of  Los  Angeles  up  to  S£uita  Barbara  eirea. 
Monitor  oil  drilling  platforms  in  Santa  Barbcura 
Channel. 

Monitor  aerosols  and  natural  pollutants  in 
vicinity  of  coal  oil  point. 

Belease  of  tracer  gas  from  ship  west  of  Santa 
Barbeira. 


The  letters  and  numbers  preceding  some  of  the  times  2u:e  codes  used 
below  for  easy  reference  to  the  various  experiments. 
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III.  Equipment  eind  Data  Acquisition 

The  R/V  ACANIA  is  126  ft.  long,  and  of  narrow  beam  emd  low  profile, 
which  makes  it  ideal  for  meteorological  work  since  disturb^ulce  of  the 
local  air  flow  is  minimal.  Sensors  are  located  at  four  levels  on  two 
masts,  one  at  the  tip  of  the  bow,  and  the  second  15  ft.  aft  of  bow.  The 
forward  mast  has  sensors  located  at  4.2  meters  and  7 meters,  those  on  the 
rearward  mast  are  at  14.7  meters  and  20.5  meters,  where  all  heights  are 
measured  edaove  the  mean  water  line.  Each  level  contains  sensors  for  de- 
tecting both  mean  and  fluctuating  parameters,  emd  on  the  top  two  levels 
the  mecin  sensors  are  0.7  meters  below  the  fluctuation  sensors. 

A sea  surface  ten^erature  sensor  was  suspended  from  a pole  which 
extended  10  ft.  beyond  the  tip  of  the  bow.  This  sensor  is  mounted  in  a 
300  gram  brass  plug  in  the  end  of  a 6 ft.  long  by  3/4  inch  piece  of  tygon 
tubing.  The  tubing  floats  and  keeps  the  sensor  on  the  surface  ( a depth 
of  approximately  1 ft.  is  averaged  because  of  bobbing  caused  by  the  ship's 
motion)  euid  also  protects  the  sensor  from  the  sea  water.  The  brass  plug 
has  a high  heat  capacity  and  smooths  fluctuations  in  temperature  that 
would  be  caused  by  the  bobbing  of  the  sensor. 

An  acoustic  sounder  for  monitoring  the  temperature  inversion  was 
mounted  on  the  ship's  fantail.  A special  enclosure  and  mounting  were  con- 
structed to  attenuate  the  rather  severe  shipboard  acoustic  noise,  which 
limits  the  usefulness  of  the  device  when  the  ship  is  underway.  The  normal 
range  of  the  sounder  is  1 kilometer,  which  is  limited  to  approximately  500 
meters  when  the  ship  is  at  full  speed. 

The  meem  sensors  at  the  four  levels  above  the  surface  are  for  wind 
speed,  temperature  emd  humidity.  Wind  speed  is  measured  with  cup  anemometers 
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which  have  a threshold  of  0.5  knots.  Humidity  sensors  are  Li  Cl  cells 
which  have  an  accuracy  of  3%.  The  temperature  sensors  cure  quartz  ther- 
mometers (including  the  sea  surface)  vrtiich  have  em  accuracy  of  0.01®C. 

The  temperature  euid  hvmidity  sensors  are  placed  in  aspirators  which  protect 
them  from  the  environment.  The  aspirators  include  radiation  shields  so  that 
the  temperature  sensors  are  protected  from  both  direct  and  reflected  radia- 
tion from  the  sun  and  also  from  heat  radiated  from  the  ship.  With  this 
system  the  precision  of  ten^erature  measurements  Is  edx^ut  0.07*C. 

Fluctuations  of  temperature  emd  wind  speed  are  detected  with  cold 
wires  emd  hot  wires  respectively.  The  cold  wires  are  2.5  y x 2 mm  platinum. 
The  hot  wires  are  60  y x 2 mm  platinm  film  on  a quartz  substrate  and  are 
operated  at  20%  overheat,  \diich  is  high  enough  to  make  them  insensitive  to 
ten^erature  fluctuations.  The  wind  speed  bridges  are  constant  ten^erature 
anemometers.  The  tenperature  bridges  are  operated  at  3 kHz  emd  very  low 
current  so  that  the  wires  £u:e  not  heated,  thus,  they  are  not  sensitive  to 
wind  speed  fluctuations.  The  hot  wires  2ure  aligned  with  their  axes  verti- 
cal so  that  they  are  sensitive  only  to  the  horizontal  conponent  of  the  wind. 

The  me2m  signals  are  averaged  for  times  that  cire  determined  by  the  con- 
ditions (usually  either  10  or  20  min.) . Acquisition,  averaging,  2uid  recording 
on  a teletype  and  magnetic  tape  are  accomplished  by  sm  NFS  developed  system, 
designated  MIDAS. 

Fluctuation  data  was  acquired  in  two  forms,  using  either  a single 
sensor,  or  paired  sensors  at  a sepeuration  of  0.3  meters.  We  used  paired 
sensors  for  tenperature  and  a single  sensor  for  wind  speed.  Paired  sensors 
result  in  spatial  filtering  of  the  signal.  When  a single  sensor  was  used, 
the  resulting  signal  was  attenuated  above  200  Hz  and  below  5 Hz  giving  tem- 
poral filtering  which  is  equivalent  to  the  spatial  filtering  with  paired 
sensors . 
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The  resultant  signals  were  recorded  on  magnetic  tape  for  later  pro- 
cessing In  the  laboratory.  The  signals  were  also  processed  to  give  the 
RMS  values  which  were  acquired  and  recorded  by  MIDAS  and  also  recorded  on 
strip  ch2u:ts.  The  strip  charts  are  used  on  shipboard  to  obtain  real  time 
e and  values. 

The  acoustic  sounder  output  Is  a strip  chart  record  of  the  height 
from  which  the  return  echo  occurs,  giving  a real  time  presentation  of  the 
Inversion  height.  This  height  Is  conpared  with  the  twice  a day  radiosonde 
results  which  identify  both  the  height  and  strength  of  the  Inversion.  Agree- 
ment to  within  20  meters  was  consistently  obtained. 

Even  though  the  R/V  ACANIA  Is  very  well  suited  to  meteorological 
measurements  Its  presence  can  disturb  the  local  air  flow  sufficiently  to 
conpromlse  acquired  data.  This  effect  can  be  reduced  sufficiently  to  be 
negligible  by  keeping  the  ship  pointed  into  the  wind.  Of  course,  this 
cannot  always  be  done  so  we  acquire  data  only  vdien  the  relative  wind  Is 
within  30®  of  the  bow.  The  second  problem  with  a ship  platform  is  motion 
due  to  roll.  Roll  can  cause  two  adverse  effects:  introducing  an  extra 
con^onent  to  the  wind  speed  fluctuations  euid  increasing  the  rotation  rate 
of  the  cup  anemometers.  Since  fluctuation  signals  are  analyzed  only  at 
frequencies  above  5 Hz  and  ship  motions  are  much  less  than  1 Hz,  there  Is 
no  observable  affect  on  the  fluctuation  signal.  Under  conditions  of  severe 
roll  (20®  roll,  4 sec.  period)  the  mean  wind  speed  measured  at  level  4 is 
elevated  by  about  1 knot  when  the  relative  wind  is  from  the  bow.  When  the 
relative  wind  is  from  90*  off  the  bow  this  effect  is  reduced  to  zero.  The 
wind  speed  results  given  in  this  report  are  not  corrected  for  this  effect. 

Conplete  descriptions  of  the  NFS  shipboard  equipment  and  the  analysis 
methods  can  be  found  in  the  references. 
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IV . Data  Reduction 


As  described  in  the  previous  section  the  data  acquired  includes: 

Profiles  of  ten^erature , humidity,  and  wind  speed,  sea  surface  temperature, 
and  records  of  the  temperature  and  wind  speed  fluctuations.  This  data  is 
reduced  to  obtain  the  following  parameters: 

Friction  velocity 

e Dissipation  rate  of  turbulent  kinetic  energy 

2 

Teii5)erature  structure  function 
D Diffusivity 

Ri  Richardson  Nxmober 

F Sensible  Heat  Flux 

n 

F Momentum  Flux 

N 

There  are  several  possible  methods  to  evaluate  these  parameters  by  utilizing 
both  the  profile  and  fluctuation  data.  Here  we  describe  only  the  method  used 
to  process  the  data  presented  in  this  report. 

All  of  the  data  collected  have  not  been  analyzed.  There  were  many  time 
periods  when  the  relative  wind  was  from  an  unfavorable  direction,  and  these 
data  were  not  2uialyzed.  Also,  there  are  periods  when  the  results  obtained 
are  obviously  in  error  and  they  have  been  discarded.  These  errors  are  due 
to  extreuieous  noise,  such  as  ship  radio  transmission  and  power  surges. 

We  use  the  profiles  to  obtain  potential  temperature  (*c) 

0 - T + 0.0098  Z , (1) 

and  the  virtual  potential  temperature 

9^  - 9 + 0.61  q T , (2) 
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emd  their  gradients  with  height.  T is  the  absolute  temperature,  Z the 
height  above  the  mean  sea  surface  in  meters,  and  q the  specific  humidity 
in  grams  of  water  vapor  per  gram  of  dry  air.  The  gradients  are  obtained 
by  fitting  6 and  6^  with  a log  profile  and  then  evaluating  the  gradients 
at  a height  of  10  deters. 

The  specific  humidity  is  found  from 

q » 6.5  X 10  ® H e3q)j^Q[A  - - C log  T]  , (3) 

where  H is  the  relative  humidity  (%) » A » 23.84,  B - 2984,  and  C = 5.03. 

In  calculating  q for  the  various  heights  we  use  a single  value  of  H 
which  is  the  average  of  the  four  measured  values.  Thus, the  dependence  of 
q on  height  comes  from  the  temperatvire  veurlation.  This  is  done  because 
the  humidity  sensors  are  not  accurate  enough  to  allow  a profile  to  be 
specified. 

Three  methods  have  been  used  to  analyze  the  fluctuation  data:  dif- 

2 2 

ference,  BMS,  ^md  spectral.  The  structure  functions  and  are 
obtained  directly  from  the  emalysls,  and  c is  found  from 


2 e 


2/3 


(4) 


In  the  difference  method  the  structure  function  is  found  by  measvuring 
the  variance  of  the  difference  in  the  variable  x at  two  points  separated 
a known  distance,  d: 

- < tx(r)  - x(r  + d)]*  > d"^^^  . (5) 

The  spectral  method  is  based  upon  the  assumption  of  "local  isotropy" 
euid  the  Kolmogorov  -5/3  slope  of  the  one-dimensional  power  spectral  density, 
^x(3c)* 
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(6) 


(J>  (k)  » 0.25  C ^ , 

X 

where  k is  the  wave  number.  Using  Taylor's  "frozen  turbulence"  hypothesis 
(k  = 2 irf/U  } we  can  find  by  performing  a fourler  spectrum  analysis  of 
a signal  In  the  frequency  domain  (f ) : 

2 2 irf  "2/3 

f - k (J)^(k)  - 0.25  ) . (7) 


Therefore , 

2/3 

Cx^-4(^)  [f^/^  (fr^(f)]  . (8) 

Here  U Is  the  meam  wind  speed  averaged  over  the  analysis  time  to  perform 
the  spectral  ainalysls.  This  expression  Is  esqpected  to  be  valid  In  the 
inertial  Subrange,  which  covers  a frequency  remge  of  approximately  0.1  to 
100  Hz. 

The  HNS  method  Is  based  upon  measuring  the  veurlance  of  the  signal 
fluctuations  between  selected  frequency  limits  t 


k 

r - <>'*rms)* 


(9) 


where  x*  Is  the  fluctuating  cosqponent  of  x.  Aie  upper  and  lower  frequency 
limits,  f^  and  f^,  are  detezmlned  by  a filter,  and  are  related  to  the 
respective  wave  numbers  by  Taylor's  hypothesis.  Using  equations  (6)  and  (9) ; 


i / Ll 

3 ' U 


2/3 


(x'  )■ 

rms 


(10) 
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The  advauntages  of  the  difference  and  BMS  methods  are  that  the  output 
voltage  ceui  be  presented  and  averaged  on  a strip  chart  for  real  time 
analysis.  The  spectral  method  Is  more  time  consuming  and  Is  normally  per- 
formed in  the  laboratory  from  tape  recorded  signals.  However,  it  has  one 
distinct  advantage  In  that  the  spectrum  can  be  viewed  and  obvious  noise 

can  be  ignored  (such  as  60  Hz  and  its  harmonics) . 

2 

For  the  data  presented  in  this  report  was  found  using  the  difference 
method,  and  e was  found  both  from  the  KMS  and  spectral  methods. 

Using  0,  0J,,  and  6^^^  the  other  peurameters  are  obtained  from  the  following 
equations : 


U*  - (k  Z 0) 


1/3 


D = k 2 Uj^  (neutral  stability) 


R.  = ^ k Z 
i T 5 


(11) 

(12) 

(13) 


Fjj  = p C (6*  U*)  , (14) 

emd 

Fm  - P (U*)*  , (15) 


where  p is  the  density  of  air  at  STP  (1.29  kg/m*),  C the  heat  capacity  of  air 

3 t 2 

(pc  = 1.31  X 10  Joule/K  m ),  g the  acceleration  of  gravity  (9.8  m/aac  ),  and 

k is  Von  constant  (0.35)  . 

Most  of  the  above  parameters  depend  on  height.  He  arrive  at  final  values 

for  all  of  the  parameters  for  an  assumed  height  of  10  meters.  To  be  specific, 

2 -2/3  -1 

0 and  0^  are  fitted  by  a log  profile,  is  fitted  by  Z , and  0 by  Z 

The  10  meter  values  for  e,  0*  and  0 * are  calculated,  then  D,  R. , F„  and  F„ 

w yw  X n M 

cure  calculated  directly. 
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V.  Results 


A.  Acoustic  Sounder  and  Radiosonde 

The  acoustic  sounder  and  radiosonde  soundings  show  the  presence 
of  a strong  marine  inversion  for  the  full  cruise.  The  height  euid  strength 
of  the  inversion  varied  with  time  and/or  location.  The  bottom  of  each 
acoustic  sounder  strip  chart  (Figures  2a  - 2d)  is  marked  with  letters 
which  correspond  to  radiosonde  graphs  (Figures  3a  - 3u)  to  show  the  times 
at  which  the  radiosondes  were  released.  Examination  of  both  sets  of  data 
shows  good  agreement,  except  for  0200  PDT  on  7/24,  for  which  the  radiosonde 
shows  the  inversion  base  approximately  150  m higher  than  the  sounder.  The 
reason  for  the  discrepancy  is  not  known  but  we  assume  that  eui  error  was  made 
in  reducing  the  radiosonde  data. 

In  Table  I we  list  the  height  of  the  base  of  the  inversion  and  the 
strength  of  the  inversion,  taken  from  the  radiosonde  data.  The  data  2u:e 
coded  with  letters  A thru  U.  Examination  of  the  graphs  shows  that  it  is 
difficult,  in  many  cases,  to  identify  accurately  the  height  emd  strength  of 
the  inversion  since  tenperature  chemges  occur  gradually.  In  general  the 
base  of  the  inversion  is  taken  to  be  the  height  at  which  the  minimum  teirp- 
erature  occurs;  the  strength  is  from  minimum  to  maximum  ten^erature  regard- 
less of  the  height  difference,  unless  there  is  a subsequent  increase  in 
temperature  identifying  a second  Inversion.  Good  exanples  of  difficult  to 
interpret  profiles  are  Figures  3b  and  3c.  There  are  several  inflections  on 
each  profile  and  interpretation  is  very  subjective. 

The  acoustic  sounder  is  most  useful  for  continually  monitoring  the 
inversion  height,  and  hence  the  mixing  depth.  Figure  2a  shows  the  begin- 
ning portion  of  the  cruise.  Including  the  initial  Investigation  o^  the 
Santa  Barb2u:a  Channel.  The  inversion  was  quite  stable  for  the  entire  time. 
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TABLE  I 


Height  and  Strength  of  the  Marine  Inversion 
t2d(en  from  Shipboard  Radiosonde  Data 

Times  are  PDT 


Code 

Date/Time 

Inversion  Height  (m) 

Inversion  Strength  (“C) 

A 

7/17/1700 

360 

14.5 

B 

7/18/0500 

0 

0.5 

240 

13. 

C 

7/19/0200 

280 

6.7 

850 

2.5 

D 

7/19/1900 

450 

0.4 

640 

9.3 

E 

7/20/0200 

0 

0.5 

680 

9.3 

F 

7/20/1900 

200 

2.1 

650 

1.8 

G 

7/21/0200 

450 

10.  + 

H 

7/21/1900 

200 

7.8 

470 

4.2 

I 

7/21/2300 

330 

0.2 

440 

7.1 

610 

4.6 

J 

7/22/0200 

240 

13.3 

K 

7/22/1900 

380 

11.9 

L 

7/23/0200 

360 

13.4 

M 

7/23/1900 

360 

9.5 

N 

7/23/2300 

520 

5.8 

700 

3.0 

0 

7/24/0200 

650 

8. 

P 

7/26/0200 

170 

8.3 

Q 

7/27/1900 

100 

12.6 

530 

1.2 

R 

7/28/0200 

110 

14.8 

S 

7/28/0500 

20 

2.8 

180 

10.5 

T 

7/28/1700 

230 

16.3 

U 

7/29/0500 

120 

13.5 

~500 

1.8 
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varying  from  200  m to  400  m.  The  exception  was  when  the  ship  rounded  Pt. 
Conception  at  0230  on  7/18.  The  presence  of  very  intense  thermal  plumes 
indicates  the  region  of  mixing  of  marine  air  and  air  that  has  passed  over 
land  and  been  warmed.  The  region  of  disturbed  air  extends  well  into  the 
channel  from  Pt.  Conception. 

From  approximately  0900  to  1300  on  7/19  an  oil  t^mker  which  was 
not  in  the  process  of  transfering  oil  was  monitored.  The  inversion  height 
was  approximately  600  m for  the  first  3 hours,  then  dropped  to  200  m during 
the  last  hour.  At  2130  of  the  S2une  day  the  monitoring  was  resumed,  but 
during  lightering  operation,  eind  continued  until  0930  on  7/20.  The  inver- 
sion height  was  600  m to  800  m.  The  lightering  operation  was  approximately 
50  kmi  from  shore  and  the  inversion  was  much  higher  than  that  found  near 
shore . 

During  the  monitoring  of  platform  Eva  the  inversion  was  below  100  m. 
The  cheinge  in  inversion  height  appears  to  have  been  a geographical  effect 
since  the  height  was  greater  both  as  the  ship  approached  and  left  the  shore- 
line. 

We  then  proceeded  to  do  a thorough  study  of  the  turbulence  in  the 
Santa  Monica  Bay,  the  study  taking  approximately  24  hours.  The  inversion 
was  at  200  m to  400  m during  all  but  the  last  5 hours,  when  it  dropped  to 
the  surface,  split  into  more  than  one  level,  and  becaune  difficult  to  inter- 
pret, (euround  1900  on  7/21)  . 

The  first  power  pleuit  SF,  release  at  0100  on  7/22  occurred  when  the 

o 

inversion  was  150  m high  and  not  sharp.  (See  Figure  2a) . Multiple  returns 
were  obtained  from  complex  thermal  structure  up  to  700  m.  At  0400  the  inver- 
sion either  dropped  to  the  surface  or  a new  inversion  rose  from  the  surface, 
reaching  an  elevation  of  250  m by  0600.  This  variation  in  the  inversion 
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height  appears  to  be  a change  with  time  rather  than  position  since  the 
ship  was  in  the  s^une  locations  at  different  times,  with  different  inver- 


sion heights  being  obtained.  The  height  emd  intensity  of  the  inversion 
remained  stable  through  the  remainder  of  the  SFg  tracer  experiment,  which 
ended  at  ~ 1200. 

The  sounder  was  not  in  operation  during  the  period  when  data  was 
being  taken  along  the  coast  from  Del  Mar  to  Long  Beach.  The  one  radio- 
sonde tedcen  during  this  period  (Figure  2e)  shows  an  inversion  height  of 
approximately  450  m. 

When  SF  was  released  from  the  power  pl2mt  the  second  time  the 
o 

inversion  height  was  approximately  600  m ( ~0200  on  7/24).  The  inversion 
remained  fairly  high  until  0800,  at  which  time  it  was  500  m,  thereupon 
it  dropped  to  200  m by  the  end  of  the  experiment  at  1100.  Again  the 
changes  that  occurred  were  teaporal  rather  than  spatial. 

The  remainder  of  the  experiments  were  performed  under  simileur  con- 
ditions. The  inversion  was  near  the  surface  (within  200  m)  and  for  much 
of  the  time  it  was  difficult  to  distinguish  an  inversion  from  thermal 
plumes . 

B . Calculated  Results 

Appendix  B presents  a table  of  the  results  that  were  conputed  from 
the  basic  data  (Tedsle  VI) . The  data  includes  dissipation  rate,  temperature 
structure  function,  turbulent  eddy  diffuslvity,  Richardson's  nuniber, 
momentum  flux,  and  heat  flux.  Figures  4a  and  4b  show  plots  of  Richardson's 
number  and  the  momentum  flux  versus  time  for  most  of  the  cruise.  The  codes 
on  the  time  axis  refer  to  the  various  e]q>eriments  listed  in  Section  II.  It 
must  be  emphasized  that  these  numbers  give  estimates  of  flux  parameters. 
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based  on  measurements  of  inertial  subrange  parameters , not  direct  measure- 
ments. 

Difficulty  Cein  be  encoiontered  in  evaluating  temperature  profiles 
to  obtain  0^,  and  9^^ . The  following  two  examples  for  temperatures  obtained 
from  the  four  levels  have  been  taken  from  Table  V and  illustrate  the  prob- 
lem. 


17.85  17.79  17.70  17.54 

17.72  17.52  17.85  17.46 


The  first  set  of  teir^ieratures  shows  a good  profile  and  it  is  easy  to  fit  6 
with  a log  profile  to  obtain  0^ . T^  for  the  second  set  (underlined)  appears 
to  be  in  error  since  it  certainly  does  match  the  gradient  shown  by  the  other 
ten^jeratures . The  reason  for  this  discrepancy  could  be  ship  influence,  or 
it  could  be  real,  this  is  unknown  but  ship  influence  is  most  likely.  How- 
ever, it  would  be  inappropriate  to  include  the  value  when  fitting  the  log 
profile,  so  the  point  was  discarded.  Level  3 temperature  data  was  not  used 
for  calculations  for  a large  percentage  of  the  data. 

Discarding  experimental  data  is  a hazardous  thing  to  do  and  raises 
questions  eibout  the  validity  of  the  final  results.  To  this  point  we  must 
recall  that  we  are  trying  to  make  flux  estimates,  and  determine  the  stability 
of  the  atmosphere  in  a land  sea  boundary  region.  The  data  we  have  taken 
indicates  that  the  land  influence  extends  well  out  to  sea  in  this  area.  We 
did  not  enco\inter  the  homogenous,  equilibrium  atmosphere  on  which  the  theory, 
and  hence  data  evaluation  methods  2ure  based.  However,  reasonable  estimates 
can  be  obtained  by  fitting  an  averaged  profile,  fmd  this  is  what  we  have 
done  by  discarding  level  3.  We  expect  a 50%  error  in  the  final  result  for 
measurements  made  this  close  to  land,  and  on  a non-stationary  platform. 
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A few  general  conclusions  can  be  drawn  for  the  cruise  as  a whole. 

Ohe  Z = 10  m value  of  Rlchudson's  number  varied  from  slightly  positive  to 
as  low  as  -1.3,  and  was  in  the  range  of  -.1  to  -.5  the  majority  of  the  time. 
Hence,  conditions  were  generally  slightly  unstable  (the  sea  surface  was 
warmer  than  the  air  within  the  first  20  meters  for  the  full  time  the  ship 
was  in  the  Southern  California  area) . Since  small  tenperature  gradients 
were  present  small  errors  in  temperature  measurements,  such  as  those  due  to 
ship  influence,  could  have  resulted  in  large  perturbations  in  the  final 
results.  This  could  account  for  some  of  the  fluctuations  evident  in  Figure 
4. 

Note  that  Table  VI  contains  a dashed  line  at  0420  on  7/26.  This  is 
due  to  the  poor  quality  of  the  data  that  was  obtained  at  siobsequent  times. 

The  waves  were  extremely  high  during  this  period;  part  of  level  1 had  been 
destroyed  by  waves  breedcing  over  the  bow,  and  all  of  level  1 was  removed. 

The  table  shows  positive  Richardson's  numbers  from  0725  to  1030,  but  the 
numbers  are  undoubtedly  too  large.  The  zeroes  indicated  for  heat  flux  are 
due  to  the  fact  that  for  Ri  > 0.2  heat  flux  will  be  zero.  The  numbers  given 
for  momentum  flux  for  this  period  are  good  estimates. 

Because  of  the  variations  in  the  data  Figure  4 shows  two  superin^sed 
representations  of  the  data.  Individual  data  points  are  shown,  which  are  fit 
by  cin  average  curve  (heavy  lines) , and  are  also  fit  by  a curve  which  shows 
fluctuations  about  the  average  (light  lines)  . Note  that  the  curve  for  fluc- 
tuations is  only  used  in  those  regions  of  the  graphs  where  fluctuations  were 
quite  rapid.  Where  the  fluctuations  can  be  easily  identified  from  the  indi- 
vidual points  no  curve  is  included  to  reduce  the  clutter  on  the  graph.  The 
average  is  an  estimate,  not  a computed  curve. 


Refering  to  Equations  11  through  15  we  see  that  all  of  the  calculated 

2 

results  depend  on  two  parameters,  6^  and  e.  {In  actuality  was  derived 

directly  from  fluctuation  data,  not  from  6^,  although  it  would  have  been 

possible  to  do  so) . Examination  of  Figure  4 illustrates  the  manner  in  which 

Ri  and  F aure  related.  A large  temperature  gradient  (unstable  case)  results 
m 

in  a large  negative  Richardson's  number  (Ri  a ^ large  momentum 

2 -2 
flux  (f^  a u^)  results  in  a small  Richardson's  number  (Ri  a u^  ) . Thus,  when 

the  flux  is  large  Ri  will  be  small  unless  compensated  by  large  temperature  or 

humidity  gradients . This  behavior  is  evident  in  the  figures . 

Not  all  of  the  fluctuation  evident  in  the  results  are  measurements 

effects.  The  most  striking  example  is  shown  from  0400  to  1100  on  7/21.  Ri 

and  F^vary  in  somewhat  of  an  oscillatory  manner,  the  variations  being  almost 

a factor  of  4.  The  values  obtained  during  this  time  are  reasonable  estimates, 

for  20  minutes  averages.  However,  a more  reasonable  estimate  of  the  fltix, 

for  example,  would  be  obtained  from  an  hourly  average.  As  has  been  pointed 

out  by  others,  any  averaging  time  shorter  than  one  hour  is  probably  too  short 

for  processing  atmospheric  turbulence  data.  We  use  the  shorter  averaging 

time  out  of  necessity  because  of  the  shipboard  platform  and  cheinges  in  course, 

position,  and  conditions.  Longer  averaging  times  could  be  used  for  those 

periods  when  conditions  remain  constant.  Because  of  the  short  experimental 

averaging  time  it  is  more  appropriate  to  use  the  average  curve  through  the 

data  shown  in  Figures  4,  than  the  individual  points  in  order  to  obtain  final 

results . 


C.  Diurnal  Variation  and  Land  Influence 

In  Figure  5 we  have  plotted  momentum  flux  versus  time  for  the  24  hour 
day,  where  data  from  all  time  periods  of  the  cruise  have  been  averaged.  The 
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numbers  that  were  used  to  obtain  this  figure  were  obtained  from  the 
average  curves  in  Figures  4,  not  from  the  individual  data  points. 

One  expects  very  little  diurnal  variation  for  open  ocezm  conditions, 
and  this  has  been  confirmed  on  previous  cruises.  In  Figure  5 we  see  a 
significant  diurnal  variation  of  the  flux. 

Data  was  taken  on  this  cruise  for  a wide  range  of  distances  from 
shore,  as  close  as  1/2  mile  and  as  far  as  30  miles  when  the  ship  was  south 
of  Catalina  Isl^uld.  In  Figure  5 we  divide  the  averaged  data  into  two 
groups,  all  data,  and  at  sea  data.  The  at  sea  data  is  that  tedcen  in  the 
neighborhood  of  Catalina  Isleuid.  The  at  sea  data  shows  a much  smaller 
diurnal  variation,  but  the  vauriation  is  large  enough  to  indicate  signifi- 
cant land  influence  in  the  area. 

In  Figure  6 we  show  the  variation  of  momentum  flux  with  wind  speed. 
The  points  on  the  graph  were  computed  directly  from  the  data  in  Table  VI 
cind  then  averaged.  When  there  was  not  much  data  availedile  at  a particular 
wind  speed,  data  from  more  than  one  speed  ware  averaged  in  order  to  improve 
the  statistics.  The  error  bars  were  obtained  from  the  square  root  of  the 
niomber  of  points  in  each  sample.  The  solid  points  are  for  data  necu:  the 
shore,  the  open  points  are  for  at  sea  (south  of  Catalina  Islemd) . The  two 
solid  lines  are  theoretical  curves  obtained  from 

2 

F^  - p [kU/£n  (Z/Z^) ] , (16) 

for  values  of  the  roughness  length  Z “0.01  emd  Z « 0.1  cm.  Previous 

o o 

measurements  over  the  open  ocean  give  values  of  Z^  near  0.1  cm  and  the 
results  presented  in  Figure  6 are  somewhat  in  agreement  with  this  value. 
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The  data  is  much  higher  than  the  prediction  for  wind  speeds  below 

5 knots.  Obviously  the  shear  produced  momentum  flux  would  be  zero  at  zero 

-2  2 

wind  speed,  but  we  measure  values  of  the  order  10  kg/msec  . The  measure- 
ments can  be  partially  explained  by  c<jnvective  flow  and  the  roll  of  the  ship. 
The  convective  flow  could  be  enhanced  by  the  ship's  presence  since  its  mean 
temperature  is  greater  than  the  water  temperature.  Quantitative  comparison 
with  the  results  are  not  feasible. 

Careful  evaluation  of  the  results  shows  that  some,  but  not  all,  of 
the  diurnal  effect  is  due  to  em  increased  wind  speed  near  land  during  those 

periods  when  high  values  of  F were  obtained.  The  data  base  reported  here  is 

m 

not  sufficiently  large  to  separate  the  effect  of  wind  speed  from  the  influence 
of  the  nearby  land  mass. 

It  is  interesting  to  note  that  these  data  imply  that  the  "Los  Angeles 
air  basin"  is  eui  area  that  extends  at  least  30  miles  to  sea. 

D.  Pollution  Model  Parameters 

Figure  Id  shows  the  course  followed  by  the  ship  to  obtain  parameters 
for  the  CARS  air  pollution  model.  Data  taken  on  the  westward  leg  on  the  lower 
part  of  the  figure  is  labeled  M2  on  Figure  4a.  The  remainder  of  the  data  is 
labeled  M3.  In  addition  the  data  labeled  Ml  was  obtained  in  the  neighborhood 
of  Catalina  Island  in  the  hope  of  finding  open  ocean  conditions  for  comparison 
purposes.  The  model  experiment  was  performed  on  7/21.  On  7/22  and  7/24  we 
performed  the  tracer  experiments  in  the  same  zurea  and  the  parameters  presented 
here  are  averages  for  all  of  the  data.  Comparing  all  of  the  data  shows  that 
the  variations  in  the  parameters  which  we  observed  were  tenoral  rather  them 
spatial . 
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We  cannot  stress  too  strongly  that  the  parameters  presented  below 
are  averages  that  were  obtained  during  a particulau:  time  of  the  yeau:,  under 
a particular  set  of  circumstances.  As  the  weather  chemges  the  parameters 
will  change.  The  following  are  our  best  estimates; 

Ri:  The  water  was  warmer  than  the  air  leading  to  slightly  unstable 

conditions.  Expected  values  au:e  in  the  range  -0.3  £ Ri  ^ -0.05. 

This  is  consistent  with  the  Richardson's  numbers  obtained  for 
other  overwater  ea^eriments  which  average  to  Ri  »=  -0.08  (open 
ocean) . 

F^;  Figures  5 and  6 can  be  used,  recalling  that  the  solid  curve  in 

_2 

Figure  5 includes  at  sea  data  and  that  values  as  high  as  6 x 10 
2 

kg/msec  were  obtained  near  Izmd  in  the  late  afternoon.  The 

diurnal/wind  speed  vauriation  is  significeuit  and  should  be  taken 

-2 

into  account.  For  sin5)lification  one  could  use  F^  * 1.5  x 10 

-2 

for  all  periods  except  1400-1800  where  a value  of  4 x 10  would 
be  appropriate.  Alternately,  Figure  6 can  be  used  to  obtain  F^  as 
a function  of  wind  speed. 

D:  The  turbulent  eddy  diffusivity  can  be  found  from  F^  using  Equations 

12  emd  15  (assuming  near  neutral  stability  and  a height  of  10  meters) . 
We  obtain: 

D - (F ^0.105)^''^  . (17) 

2 

Near  land  this  gives  the  range  0.3  £D  _<  0.8  m /sec.  This  is  much 
lower  than  overland  values. 

e:  The  dissipation  rate  can  also  be  found  easily  frcn  F^.  Using  Equa- 

tions 11  and  15,  for  a height  of  10  meters,  we  obtain: 
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(18) 


e = (F  /2.97)^/^  . 
rn 

F : The  sensible  heat  flux  varies  quite  widely  since  it  depends  both 
H 

on  the  turbulence  and  on  the  temperature  profile.  Of  course  F 

H 

can  be  either  positive  or  negative  depending  on  the  direction  of 

the  gradient.  We  obtained  values  that  rsmged  from  -50  to  +20 

2 2 
watts/m  . A reasonable  average  would  be  -4  watt/m  for 

-2 

F^  ~ 1.5  X 10  and  Ri  •“  -.3.  For  the  same  stability  a momentum 

-2  2 
flux  of  6 X 10  will  give  a heat  flux  of  ^ -20  watt/m  . For  a 

_2 

much  larger  thermal  gradient  (Ri  ^ -3)  a momentum  flux  of  ■■  1 x 10 
will  produce  the  same  heat  flux.  Further  extrapolations  can  be 
made  but  it  is  best  to  refer  to  Table  VI. 

E.  Stability  Corrections 

As  has  been  described  above,  the  data  was  reduced  assuming  neutral 
stability.  Conditions  were  slightly  unstable,  the  average  Richardson's  nxxn- 
ber  being  -0.3.  This  value  of  Richardson's  number  leads  to  a stability  cor- 
rection of  20  to  30%  for  the  calculated  results,  depending  on  the  parameter. 
This  correction  is  within  the  errors  associated  with  the  measurements  and 
has,  hence,  not  been  applied. 
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APPENDIX  A:  BASIC  DATA 

In  this  appendix  we  present  data  on  a)  the  ship's  speed  and  course 
and  cruise  charts,  b)  wind  speed  and  direction,  c)  acoustic  sounder  strip 
charts  and  interpretation,  d)  radiosonde  results,  e)  sea  surface  tempera- 
ture and  air  temperature  and  humidity. 

A , Ship  Course 

Table  II  lists  changes  in  the  ship's  heading  emd  speed,  the  informa- 
tion being  obtained  from  the  ship's  log  emd  the  NPS  scientific  log.  The 
ship's  normal  cruising  speed  is  approximately  9 1/4  knots  fuid  is  designated 
as  full  ahead  in  the  table.  No  correction  is  made  for  currents  or  winds. 
Charts  of  the  ship's  position  are  shown  in  Figures  la  - 1 j . Small  course 
changes  which  are  occasionally  necesscury  for  maneuvering  purposes  are  not 
shown  in  either  the  table  or  the  charts. 
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TABLE  II 


IIVIIII 


i/i%m 


7119/77 


Changes  in  Ship's  Speed  and/or  Course 


0815  Departed 

0850  c/c  210“ 

0924  c/c  168“ 

1126  c/c  150“ 

1400  c/c  148“ 

1500  c/c  143“ 

1655  c/c  144“ 

0100  c/c  126“ 

0220  c/c  079“ 

0335  Downwind  leg  of  platform  Helen  run 

0345  Start  360“  pass  of  platform 

0355  Underway  to  Holly  c 100“ 

0615  Stop  for  true  wind,  maneuvering  near  Holly 

0654  c/c  105*  Depart  Holly 

0815  Stop  for  true  wind,  1.5  mi  from  platform  C 

0822  c 090“ , 1/2  ahead  near  platforms 

0835  c/c  145“  into  clean  air 

0855  c/c  070“,  1/2  ahead  run  by  platform 

0915  c 133“,  full  ahead 

1135  c/c  115“ 

1230  c/c  110“ 

1355  c/c  104“ 

1550  Drifting 

1630  Head  into  wind,  1 engine  dead  slow 

1845  c/c  060“ 

1930  c/c  270“ , port  engine  slow  ahead 

2242  c/c  080“ 

2350  c/c  270“,  starboard  engine  slow  ahead 


0215  Underway  to  beach 

0255  Drifting 

0345  Underway  to  Catalina 

0500  c/c  148“ 

0750  c/c  205“ 

0900  maneuvering  near  Edinburgh 

0945  c/c  345“,  full  ahead 

1120  Downwind  Edinburgh,  160“  at  2 knots 

1140  Finished  downwind  leg 

1149  Upwind,  160*  at  2 knots 

1205  Finished  tanker  data,  underway  to  Isthmus 

1425  In  fisherman's  cove 

1455  Underway  easterly 

1545  Heading  into  wind 

1730  Underway  to  tanker 

1800  c 133*,  full  ahead 

1935  c/c  165“ 

2120  Drift  neighborhood  of  tanker 
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7/21/77 


7/22/77 


0820 

c 348°,  full  ahead 

1015 

c/c  345° 

1445 

Underway  to  platform  Eva 

1710 

Arrive  at  Eva 

1800 

c 210°,  half  ahead 

1815 

c/c  270° 

2200 

Stop  for  true  wind 

2205 

c 030°,  full  ahead 

2300 

Stop  for  true  wind 

2325 

c 000°,  full  ahead 

0010 

c/c  345° 

0100 

Stopped  for  true  wind 

0250 

Stopped 

0315 

Underway,  c 180° 

0345 

Stopped 

0410 

Underway,  c 177* 

0445 

Stopped 

0510 

Underway,  c 175* 

0540 

Stopped 

0606 

c 180°,  full  ahead 

0635 

Stopped 

0640 

c 240°,  dead  slow 

0700 

Stopped 

0810 

c 090°,  2/3  ahead 

0900 

c/c  000° 

1110 

Stopped 

1245 

Underway,  c 265°,  1/2  ahead 

1410 

c/c  182°,  full  ahead 

1440 

At  position,  head  into  wind 

1515 

Underway  to  new  position,  full 

ahead 

1555 

At  new  position,  head  into  wind 

1620 

Underway  to  new  position,  full 

ahead 

1700 

At  position,  head  into  wind 

1725 

Underway  to  new  position,  full 

ahead 

1800 

At  position 

1823 

Underway,  c 078° 

1925 

At  position,  dead  slow  into  wind,  c 270° 

1947 

Underway,  c 000°,  full  ahead 

2011 

At  position,  c 235°,  dead  slow 

into  wind 

2038 

Full  ahead,  c 010° 

2110 

At  position 

2135 

Full  ahead,  c 053° 

2212 

Stop  for  true  wind 

2217 

Resume  course  and  speed 

2305 

Stop  for  true  wind 

2310 

c 000°,  full  ahead 

2333 

c/c  180°,  1/2  ahead 

2350 

Stopped 

0001 

c 350° 

0050 

c 009° 

0055 

c 026° 

0115 

c/c  168° 

0150 

c/c  208° 

0200 

c 200° 

A-3 


mini 

(cont'd) 


l/2im 


0211  c/c  336* 

0222  Underway,  c 336* 

0258  c/c  023* 

0330  c/c  012* 

0350  c/c  270* 

0415  Stop  for  true  wind 

0420  Underway,  c 216*,  full  ahead 

0450  c/c  200* 

0510  Stop  for  true  wind 

0530  Underway,  c 135* 

0615  Stop  for  true  wind 

0634  c 000*,  full  ahead 

0710  Stop  for  true  wind 

0756  Stop  for  true  wind 

0759  c 000*,  1/2  ahead 

0805  c/c  150*,  full  ahead 

0818  c/c  170* 

0831  2/3  ahead 

0835  c/c  090*,  1/2  ahead 

0845  Full  ahead 

0850  c/c  180* 


0910  c/c  270* 

0945  Stop  for  true  wind 

0948  c 260*,  full  ahead 

1005  Stop  for  true  wind 

1040  Stop  for  true  wind 

1044  c 122*,  full  ahead 

1130  Stop  for  true  wind 

1150  c 085* 

1305  Rendevous  at  Redondo  Beach 

1315  Underway  to  Catalina,  c 215* 

1405  c/c  155* 

1700  c/c  180* 

1800  Drifting  near  Edinburgh 

2100  Underway,  c 100*,  full  ahead 


0100  c 095* 

0200  c/c  085* 

0235  Stop  for  true  wind 

0245  Underway,  c 325*,  full  ahead 

0500  c 300* 

0530  d/c  304* 

0625  c/c  292* 

0845  c/c  240* 

0915  c/c  305* 

1130  Reduce  speed  to  slow  ahead 

1150  Drifting 

1420  Ship  stationary  pointing  into  wind 

1450  Underway,  c 215*,  1/2  speed  into  wind 

1510  Full  ahead  into  wind 

1530  c 035*,  full  ahead 

1550  c 000* 

1630  Stop  for  true  wind 

1650  No  motion,  point  into  wind 


A-4 


1720 

c 265°,  1/2  ahead  into  wind 

1745 

c 265°,  full  ahead  into  wind 

1810 

Stop  for  true  wind 

1812 

c 025°,  full  ahead 

1820 

c/c  100° 

1932 

c/c  090° 

1943 

Stopped 

2315 

Underway,  c 316°,  1/2  ahead 

l/lMii 

0100 

Drifting 

0110 

Underway,  c 092°,  1/2  ahead 

0150 

Dead  in  water 

0210 

Move  closer  to  beach 

0240 

Underway,  c 245°,  1/2  ahead 

0300 

c 245°,  full  ahead 

0450 

Stop  for  true  wind 

0455 

Underway,  c 050°,  full  ahead 

0530 

c/c  038° 

0548 

Stop  for  true  wind 

0600 

c 150°,  full  ahead 

0707 

Stop  for  true  wind 

0745 

Underway,  c 150°,  full  ahead 

0855 

Stop  for  true  wind 

0905 

c 055°,  full  ahead 

0946 

Stop  for  true  wind 

0951 

Underway,  c 330° 

1112 

Stop  for  true  wind 

1130 

Underway,  c 265°,  full  ahead 

1140 

c/c  220° 

1150 

c/c  100° 

1245 

Drifting 

1330 

Underway  to  LA  Seabuoy 

1800 

Maneuvering  in  LA  harbor 

nism 

1900 

Depart  pier  LA  harbor 

1938 

Clear  harbor  entrance,  c 198 

2027 

Stop  engines 

mtm 

0100 

Drifting 

0345 

Underway  to  Pt  Fennin 

0445 

Arrive  at  station 

0530 

Underway,  c 270°,  full  ahead 

0630 

c/c  300° 

0930 

Stop  for  true  wind 

1200 

c 298°,  speed  9 knots 

1210 

c/c  285° 

1310 

Stop  for  true  wind 

1330 

c/c  278° 

1425 

Reduce  speed  to  1/2  ahead 

1430 

Resume  speed 

1501 

Stop  for  true  wind 

1530 

Reduce  speed  to  8 knots 

1645 

Dead  slow 

1650 

Resume  full  ahead 

1700 

Speed  8 knots 

1735 

c/c  080°,  1/2  ahead 

1800 

c 075°,  slow  ahead 

2312 

1/2  ahead  for  oil  platforms 

2400 

Nearing  platforms 
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i/nni 

0448 

Complete  working  platforms 

Underway,  c 185*,  full  ahead 

0535 

Drifting 

1145 

Underway,  c 122",  full  ahead 

1345 

Arrive  channel  island  harbor 

1630 

Underway  for  Coal  Oil  Point 

1800 

c 292",  full  ahead 

2020 

Stop  for  true  wind,  near  Holly 

2235 

Leave  Coal  Oil  Point  area,  c 180",  slow  ahead 

nizm 

0100 

Drifting 

0420 

Underway  to  station,  c 245" 

0500 

Drifting  on  station 

0719 

Underway,  c 292" 

1055 

c/c  303" 

1225 

c/c  330" 

1255 

Heavy  seas,  3/4  ahead 

1435 

Heavy  seas,  1/2  ahead 

7129/11 

0245 

c/c  320" 

0340 

c/c  325",  speed  5 knots 

0600 

c/c  349" 

1200 

At  Cypress  Point 

A-6 


•uino 


Dume 


0200 


A-13 


1 1 8“  30' 


Santa 


B.  Wind  Speed  and  Direction 

Table  III  lists  ship's  heading  emd  speed,  relative  wind  direction 
and  speed , and  true  wind  direction  amd  speed . All  speeds  are  in  knots , 
and  directions  are  in  degrees  measured  clockwise  from  true  North  for  the 
ship  heading  eind  true  wind,  and  clockwise  from  the  ship's  bow  for  the 
relative  wind.  The  ship  speed  was  obtained  from  the  bridge  and  is  based 
on  the  speed  of  the  ship's  screws.  Relative  wind  was  measured  both  from 
the  ship's  anemometer  and  from  the  cups  on  level  4 of  the  NFS  scientific 
equipment.  True  wind  was  calculated  from  the  other  data  and  is  no  more 
accurate  than  ^ 1 knot  and  ^ 20®,  especially  at  low  wind  speeds. 
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TABLE  III 


Relative  and  True  Wind 


Ship 

Relative 

Wind 

True  Wind 

Time 

Heading 

Speed 

Direction 

Speed 

Direction 

Speed 

7/17 

1000 

168 

10 

013 

12.5 

220 

4 

1100 

168 

10 

010 

9.5 

279 

2 

1200 

150 

10 

010 

6 

316 

4 

1300 

150 

10 

030 

5 

306 

6 

1400 

148 

10 

020 

5 

310 

6 

1500 

143 

9 

038 

3 

307 

7 

1600 

143 

9 

087 

2 

310 

9 

1700 

144 

9 

092 

3.5 

303 

10 

1800 

144 

9 

082 

4.5 

296 

9.5 

1900 

144 

9 

090 

4 

300 

9.8 

2000 

•144 

9 

006 

2 

322 

7.0 

2100 

144 

9 

349 

1.5 

326 

7.5 

2200 

144 

9 

on 

3 

319 

6 

7/18 

0000 

144 

9 

300 

5.5 

001 

8 

0300 

079 

9 

001 

4 

258 

5 

0330 

- 

295 

4 

0400 

100 

9 

012 

8.5 

211 

2 

0500 

102 

9 

028 

9 

206 

4 

0600 

075 

9 

000 

16.5 

075 

7.5 

0630 

- 

075 

9 

0800 

097 

9 

023 

13.5 

154 

6 

0830 

095 

2.5 

020 

7 

125 

5 
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7/18 


0900 

070 

3 

020 

6.5 

106 

4 

1000 

133 

9 

032 

6.2 

272 

5 

1100 

133 

9.5 

047 

3 

297 

8 

1200 

113 

9.5 

047 

4 

270 

7 

1300 

111 

9.5 

019 

6 

264 

4 

1400 

104 

9.5 

031 

2 

276 

8 

1600 

- 

265 

7 

1700 

277 

2 

063 

10 

351 

9 

1800 

275 

2 

075 

8.5 

004 

8 

1900 

073 

9 

079 

3 

234 

9 

2000 

264 

2 

345 

8 

244 

6 

2100 

265 

2 

340 

3.5 

222 

2 

2120 

265 

2 

000 

4 

265 

2 

2140 

265 

2 

000 

5 

265 

3 

2220 

260 

2 

315 

4 

186 

3 

7/19 

0000 

265 

2 

010 

4.5 

283 

3 

0020 

270 

2 

020 

4 

308 

2 

0040 

270 

2 

018 

3 

317 

1 

0100 

270 

2 

018 

3 

317 

1 

0120 

265 

2 

060 

2 

025 

2 

0140 

275 

2 

012 

2 

Oil 

0.5 

0200 

270 

2 

350 

3 

241 

1 

0220 

090 

9 

270 

8 

312 

1 

0300 

194 

3 
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7/19 

0400 

180 

9 

000 

11 

180 

2 

0440 

180 

9 

355 

12 

161 

3 

0500 

150 

9 

000 

10 

150 

1 

0520 

148 

9 

004 

11.5 

166 

3 

0600 

140 

9 

350 

9 

045 

2 

0700 

140 

9 

355 

9 

048 

1 

0900 

- 

200 

4 

1000 

353 

2 

226 

3 

201 

5 

1100 

175 

9 

357 

11.5 

161 

3 

1140 

155 

2 

017 

3.5 

192 

2 

1200 

168 

2 

000 

5 

168 

3 

1330 

305 

9 

000 

0 

125 

9 

1430 

254 

0 

340 

9 

234 

9 

1540 

- 

328 

4 

1550 

320 

2 

000 

7 

320 

5 

1610 

320 

2 

340 

8 

294 

6 

1630 

290 

2 

355 

8 

283 

6 

1650 

285 

2 

350 

9.5 

272 

7.5 

1710 

275 

2 

000 

9 

275 

7.0 

1800 

133 

9 

058 

2 

301 

8 

1900 

133 

9 

327 

4.5 

338 

6 

2000 

167 

9 

053 

10 

277 

8.5 

2100 

147 

9 

028 

6 

290 

5 

2120 

- 

005 

3 

2140 

350 

0 

328 

3 

318 

3 

2300 

315 

'^O 
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7/20 


0000 

- 

183 

1 

0100 

- 

067 

4 

0200 

- 

290 

4 

0300 

- 

090 

6 

0400 

- 

280 

7 

0500 

- 

280 

6 

0555 

- 

270 

5 

0640 

295 

0 

000 

6 

295 

6 

0700 

- 

285 

7 

0830 

345 

9 

340 

7 

210 

3 

0850 

350 

9 

340 

9 

250 

3 

0920 

350 

9 

335 

9 

248 

4 

1000 

350 

9 

315 

10 

245 

7 

1020 

345 

9 

340 

7.5 

218 

3 

1100 

343 

9 

040 

6.5 

117 

6 

1200 

340 

9 

345 

5 

177 

4 

1220 

340 

9 

340 

6 

193 

4 

1300 

332 

9 

338 

7.5 

206 

3 

1400 

- 

280 

4 

1500 

152 

9 

068 

9.5 

274 

10 

1600 

102 

9 

181 

6.5 

282 

15 

1700 

094 

9 

146 

5.5 

261 

14 

1710 

- 

260 

12 

1720 

00 

0 

068 

4 

255 

4 

1740 

350 

0 

285 

11.5 

275 

11 
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7/20 

1800 

210 

5 

342 

19 

186 

14 

1820 

270 

5 

000 

17 

270 

12 

1840 

275 

5 

359 

16 

274 

11 

1900 

275 

5 

000 

17 

275 

12 

1920 

250 

5 

000 

14 

250 

14 

1940 

270 

5 

000 

14 

270 

14 

2000 

270 

5 

000 

16 

270 

11 

2020 

270 

5 

000 

15 

270 

10 

2045 

- 

280 

7 

2050 

270 

5 

000 

11 

270 

6 

2100 

270 

5 

000 

14 

270 

9 

2120 

270 

5 

000 

12 

270 

7 

2200 

- 

280 

4 

2300 

- 

302 

4 

2325 

000 

9 

000 

6 

180 

3 

7/21 

0000 

000 

9 

330 

9 

255 

5 

0040 

345 

9 

333 

10.5 

259 

5 

0100 

- 

305 

3 

0120 

270 

9 

350 

8 

141 

2 

0140 

270 

9 

348 

6 

112 

3 

0200 

270 

9 

350  ■ 

5 

102 

4 

0215 

- 

120 

5 

0250 

- 

100 

5 

0345 

100 

2 

000 

7 

100 

5 

0350 

095 

0 

350 

6 

085 

6 

0447 

- 

142 

5 

Pi-22 


0500 

- 

125 

4 

0545 

- 

160 

3 

0635 

- 

calm 

0830 

095 

7 

350 

12 

072 

5 

0850 

095 

7 

360 

13 

095 

6 

0900 

360 

7 

030 

8 

091 

4 

0930 

360 

9 

010 

8 

129 

2 

1000 

360 

9 

020 

7 

135 

3 

1030 

330 

9 

000 

9 

*^0 

1050 

330 

9 

350 

9 

235 

1 

1115 

- 

270 

7 

1220 

- 

270 

10 

1245 

265 

5 

000 

16 

265 

11 

1300 

265 

5 

355 

19 

258 

14 

1400 

'260 

5 

018 

18 

285 

13 

1445 

- 

280 

15 

1600 

- 

260 

14 

1700 

- 

270 

11 

1800 

- 

265 

11 

1925 

- 

250 

1 

8 

2010 

- 

225 

9 

2115 

- 

220 

2 

2215 

- 

120 

4 

2300 

- 

100 

3 

2320 

000 

9 

015 

10 

076 

3 

2350 

_ 

050 

2 
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7/22 


0004 

- 

040 

2 

0030 

- 

010 

1 

0045 

- 

025 

2 

0100 

- 

020 

2 

0130 

170 

9 

340 

10 

087 

3 

0211 
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9 
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9 
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C . Acoustic  Sounder 

Reproductions  of  the  acoustic  sounder  strip  charts  are  shown  in 
Figures  2a  - 2d.  Each  strip  in  the  photographs  shows  a 24  hour  period 
and  the  vertical  span  is  1000  meters.  Return  echos  are  readily  apparent 
and  in  most  cases  one  strong  return  identifies  the  presence  of  a temperature 
inversion.  The  height  of  the  base  of  the  inversion  is  identified  on  the 
charts  from  the  lower  most  portion  of  the  return  echo  signal  (dark  horizon- 
tal beind  on  chart)  . The  vertical  d^k  areas  on  the  charts  are  caused  by  the 
increased  noise  when  the  ship  is  in  motion.  Thermal  plumes  are  evident  on 
many  of  the  charts  and  can  be  identified  by  a return  at  the  surface  cuid 
extending  upward,  at  times  to  heights  of  nearly  300  meters.  The  times  at 
which  radiosondes  were  launched  are  shown  at  the  bottom  of  each  chart  by  a 
letter,  which  corresponds  to  the  letter  designation  on  the  radiosonde  graphs. 
Experiment  codes  are  also  shown. 

Table  IV  is  a compilation  of  data  tedcen  from  the  acoustic  sounder  charts. 
We  present  this  table  because  it  is  much  more  difficult  to  identify  weak 
returns  on  the  reproductions  of  the  charts  than  on  the  originals,  and  at 
times  as  m^my  as  5 or  6 distinct  returns  ceui  be  identified.  The  table  lists 
the  height  of  top  of  thermal  plumes,  when  present,  and  the  height  of  the 
base  of  non-surface  return  regions,  all  heights  being  in  meters.  There  eure 
also  code  letters  associated  with  the  listed  heights.  The  letters  "W"  euid  "S" 
refer  to  weak  and  strong  return  echos,  respectively,  when  they  precede  a 
number.  A "W"  following  the  listed  height  for  a return  height  Indicates  that 
echos  aure  obtained  over  a wide  band.  "Dark"  means  that  there  is  too  much 
background  noise,  giving  a dark  trace,  to  allow  determination  of  the  presence 
of  a return. 
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Figure  2c 


TABLE  IV 


SUMMARY  OF  ACOUSTIC  SOUNDER  RESULTS 


DATE/ 

THERMAL 

RETURN 

HEIGHTS 

TIME 

PLUMES 

1 

2 

3 

4 

COMMENTS 

7/16/77 

1120 

S290 

no  thermal  plumes  untl  1400 

1140 

S360 

1200 

S370 

1220 

S340 

460 

divides 

1240 

S380 

420 

1300 

S400 

460 

1320 

W300 

380 

480 

return  at  300m 
appears  - 2nd  return 

1340 

W300 

W400 

480 

getting  weak  - 3rd 
return  rising,  weak 

1400 

260 

S360 

7 

480 

return  182  merge 

1420 

200 

380 

? 

510 

can't  tell  where  #2  begins 

1440 

W160 

360 

? 

520 

1500 

220 

300 

7 

W520 

seems  to  be  4 returns  but 
can't  tell  where  it  starts 

1520 

160 

300 

7 

540 

1540 

200 

260 

7 

520 

1600 

160 

S320 

420 

return  384  merge  together 

1620 

160 

S300 

400 

500 

1640 

W130 

280 

merged 

W480 

W580 

1700 

250 

320 

500 

600 

1720 

240 

300 

merged 

520 

580 

1740 

260 

— 290 

-360 

VW540 

plumes  merge  with  1st 
return 

1800 

240 

320 

- 400 

530 

return  #2  merges  with  1 

around  400  m 


A-37 


1820 

290 

360 

dark 

560 

620 

1840 

120 

180 

300 

— 360 

560 

1900 

SlOO 

160 

320 

VW580-600 

1920 

SlOO 

180 

300 

340 

W580 

may  be  another  return  at 

460  m 

1940 

140 

—200 

260 

#3  merges  with  2 

2000 

140 

260 

merged 

340 

2020 

VW  100 

140 

280 

320 

VW420 

inversion  rising  from 

surface 

2040 

It 

140 

260 

320 

If 

return  3 seems  to  have 
split 

2100 

•^120 

240 

300 

340 

400 

2120 

150 

200 

410 

may  have  5 return  heights 

2140 

240 

combined 

VW450 

3 is  combined  with  2 and 
another  one 

2200 

160 

290 

340 

- 420-460 

return  4 is  hazy  - not  clear 

2220 

160 

290 

— 340 

440 

2240 

260 

— 310  -330 

merged 

W460 

no  plumes 

2300 

VW200 

290 

-300  - 

400 

460 

2320 

180 

260 

return  2,3,4  merged 

2340 

220 

260 

2400 

180 

S270 

360 

7/17/77 

0020 

W200 

280 

400 

Inversion  rising  from 

surfaces 

0040 

230 

S280 

380 

no  plumes 

0100 

280 

400 
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0120 

W80 

320 

420 

can't  discern  return 

1 8 3 

0140 

WlOO 

280 

310 

-'410C?) 

0200 

260 

280  merged 

0220 

WlOO 

290 

460 

there  are  now  4 returns 
all  merging 

0240 

W120 

260 

not  clear 

0300 

W160 

240 

no  5 returns 

0320 

W140 

220 

300 

0340 

240 

280  320 

~400 

0400 

220 

460 

0420 

240 

merged 

330 

inversion  from  surface 

0440 

80 

220 

280 

320 

#5  hazy 

0500 

120 

220 

270 

310 

485  merged  together 

0520 

140 

250 

310 

0540 

160 

260 

340 

0600 

60 

260 

310  420 

0620 

80 

280 

merged 

can't  tell  #384 

0640 

160 

280 

320 

0700 

200 

0720 

260 

300 

no  plumes 

0740 

100 

S200 

340 

0800 

160 

210 

. 390  350 

0820 

200 

300 

350 

0840 

-'210 

300 

340 

0900 

240 

getting  too  dark  to 
see  much 

0920 

320 

0940 

1000 
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1020 


280 


340 


1040 

1100 

300 

1120 

200V 

1140 

W80 

200  ■ 

1200 

100 

180N 

240 

1220 

100 

160W 

124C 

120 

180 

1300 

120 

220 

1320 

120 

210 

1340 

130 

220 

1400 

120 

220 

1420 

140 

210 

1440 

215 

o 

260 

1500 

130 

240 

300 

1520 

90 

240 

300 

1540 

110 

240 

305 

1600 

120 

210 

300 

1620 

140 

270 

320 

1640 

no 

280 

1700 

140 

280 

1720 

no 

300 

1740 

130 

300 

310 

1800 

120 

260 

360 

1820 

140 

280 

380 

1840 

160 

300 

1900 

140 

320 

may  have  2nd  return 


very  faint  second  return 


have  2nd  return  but  can't 
judge  - too  dark 
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1920 


140 


300 


1940 

140 

280 

2000 

150 

300 

2020 

140 

300 

2040 

260 

280 

2100 

160 

320 

2120 

120 

310 

2140 

200 

300 

2200 

160 

300 

2220 

150 

300 

2240 

130 

280 

2300 

180 

280 

2320 

120 

300 

2340 

200 

290 

2400 

160 

280 

7/18/77 

0020 

140 

260tf 

0040 

120 

260 

0100 

120 

255 

0120 

125 

145 

0140  can't  tell  - too  dark 

0200 
0220 

0240  160(7)  220W 

0300 


return  1 descending  to 
surface 

may  be  2 or  3 returns 


inversion  - going  to 
surface  - plumes  become 
darker,  rising  up  blotch 
so  can't  tell  which  is 
which 

may  have  2 or  3 returns 
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0320 

140 

160W 

0340 

140 

180 

400 

0400 

100 

180 

360 

inversion  rising  from 
surface 

0420 

120 

140 

240 

plumes  rising  to  combine 
with  inversion  at  120-140m 

0440 

80 

HOW 

240 

too  dark  to  see  3rd 

0500 

160 

220 

0520 

120 

160 

0540 

80 

100 

250 

0600 

140 

220 

0620 

120 

200 

0640 

180 

200 

310 

0700 

too  dark  to  see  anything 

0720 

0740 

0800 

0820 

120 

160 

0840 

130 

200 

0900 

100 

240 

340 

0920 

90 

210 

340 

0940 

100 

240 

2nd  return  may  have 
combined  with  1st 

1000 

100 

220 

1020 

120 

200 

1040 

140 

160W 

1100 

130 

180 

1120 

150 

200 

1140 

160 

230 

1200 

150 

260 

A-42 

I 


2120 


WlOO 


420 


2140 

WlOO 

420 

2200 

120 

380 

2220 

120 

380 

2240 

100 

380 

2300 

110 

340 

2320 

320 

2340 

140 

340 

2400 

160 

280 

7/19/77 

0020 

140 

330 

0040 

300 

0100 

320 

• 

0120 

W150 

270W 

0140 

W150 

190W. 

360 

• 

0200 

120 

240W 

480 

S20 

0220 

160 

240W 

OAMC 

0240 

100 

DARK 

0300 

120 

180 

380 

0320 

130 

140 

340W 

460 

540 

0340 

140 

-'US 

370W 

0400 

100 

360 

DARK 

0420 

TOO 

DARK 

0440 

0500 

0520 


dividing  into  3 returns 

pluses  are  very  faint 

can't  tell  where  3 returns 
start 


inversion  descending  to 
surface 

3 I 4 merged  with  2 

too  dark  to  see  until  0850 


0540 
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0600 

0620 

0640 

0700 

0720 

0740 

0800 

0820 

0820 

0840 


0900 

220 

S580 

only  1 inversion  now 

0920 

220 

590 

0940 

240 

620 

1000 

240 

370 

580 

another  return  appeared 
at  370  m 

1020 

160 

300 

620 

1040 

DARK 

1100 

DARK 

1120 

140 

W240 

610 

1140 

210 

at  surface 

600 

1st  return  at  surface 

1200 

180 

tf 

600 

1220 

180 

tl 

580 

1240 

160 

disappeared 

420 

#2  dropping 

1300 

200 

260 

almost  at  surface 

1320 

200 

220 

1340 

160 

300 

1400 

140 

340 

A-45 


inversion  very  wide  not 
sure  if  it's  part  of 
plume  or  inversion  (200-660) 


1420 

170 

1440 

190 

1500 

180 

1520 

160 

1540 

180 

1600 

190 

1620 

210 

1640 

220 

1700 

220 

1720 

220 

1740 

240 

1800 

220 

1820 

160 

1840 

160 

1900 

160 

1920 

160 

1940 

180 

2000 

220 

2020 

240 

2040 

200 

2100 

160 

2120 

180 

2140 

200 

2200 

180 

2220 

140 

2240 

200 

2300 

130 

200W 


200N 

TOO  DARK 


If 

500 

500 

470 

60« 

500 

500 

400 

568 

420 

510 

410 

MIK 

460 

SIO 

- 

540 

DARK 

500 

500 

511 

540 

560 

sto 

600 

590 

600 

580 

650 

2nd  return  appeared 
#2  merfed  with  #1 
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2320 

150 

650 

2340 

2400 

200 

640 

7/20/77 

0000 

140 

620 

0020 

140 

635 

0040 

180 

630 

0100 

180 

620 

0120 

140 

640 

0140 

140 

640 

0200 

120 

690 

0220 

200 

720 

0240 

320 

750 

0300 

380 

760 

0320 

380 

790 

0340 

400 

760 

0400 

380 

740 

0420 

360 

720 

0440 

360 

740 

0500 

300 

720 

0520 

320 

700 

0540 

280 

280 

690 

0600 

can't  tell  where 
plume  ends  ft 
inversion  begins 

709 

0620 

180 

120 

690 

0640 

7 

160 

620 

0700 

7 

160 

615 

another  return  appears  at 
280  m 

very  faint  - may  have 
merged  with  plumes 
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0720 


160 


220 


615 


0740 

180 

230 

600 

' 

0800 

200 
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D . Radiosonde 


Figures  3a  - 3u  show  the  radiosonde  results  up  to  elevation  of  1000 
meters.  Results  eure  not  shewn  above  this  height  because  we  are  prinuurily 
interested  in  the  marine  Inversion.  Ihese  figures  were  derived  directly 
from  the  radiosonde  readout  strip  chart  presentation  so  that  fine  detail 
can  be  observed.  A good  exanqple  is  shown  in  Figure  3b  ^ere  the  results 
show  a surface  Inversion « an  elevated  inversion  with  a 230  m base,  and  10 
changes  in  slope  above  the  inversion.  Detail  of  this  type  is  not  avail- 
able in  the  normal  presentation  of  radiosonde  results  where  only  major 
features  are  plotted  on  a much  larger  height  scale. 
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E.  Tenperature  and  Humidity  Data 

Table  V lists  sea  surface  ten^erature , ad.r  ten^rature  at  the  four 
shipboard  levels  and  the  average  relative  humidity  from  the  four  levels. 
Due  to  inherent  inaccuracies  in  the  humidity  sensors  it  would  be  inapprop- 
riate to  list  the  four  values  obtained  from  the  levels . All  tenperatures 
are  in  “C.  The  only  data  presented  in  Table  IV  is  for  those  time  periods 
iidien  it  was  possible  to  obtain  good  profiles  and  fluctuation  data  (those 
times  when  the  relative  wind  was  from  a favorable  direction) . Data 
entries  are  missing  for  some  times  due  to  malfunctions  of  the  teletype  or, 
in  the  case  of  the  sea  surface  and  level  1 for  the  latter  part  of  the 
cruise,  due  to  system  damage  by  high  waves. 
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TABLE  V 


Sea  Surface  Teoperature  and  Atmospheric  Tenperature 
and  Relative  Humidity 


Time 

^1 

^2 

^3 

•^4 

H 

7/18 

2025 

16.41 

16.37 

16.34 

16.25 

90 

2120 

17.11 

16.45 

16.42 

16.36 

16.27 

90 

2145 

16.49 

16.49 

16.47 

16.42 

16.33 

90 

2200 

16.56 

16.53 

16.49 

16.40 

89 

2220 

16.42 

16.51 

16.47 

16.42 

16.34 

90 

2245 

16.47 

16.42 

16.33 

90 

7/19 


0000 

19.11 

16.53 

16.44 

16.36 

90 

0020 

18.54 

16.34 

16.47 

16.28 

90 

0100 

16.13 

16.08 

16.06 

16.01 

92 

0140 

16.70 

15.80 

15.89 

15.94 

15.90 

93 

0245 

16.06 

16.03 

15.98 

15.95 

15.92 

100 

0420 

17.04 

16.58 

16.56 

16.51 

16.43 

90 

0440 

17.4 

16.7 

16.67 

16.62 

16.53 

90 

0500 

17.70 

16.71 

16.63 

16.33 

89 

0520 

17.7 

16.8 

16.73 

16.67 

16.58 

87 

0600 

16.81 

16.92 

16.90 

16.89 

16.80 

85 

0620 

19.20 

17.31 

17.08 

84 

0640 

19.00 

17.51 

17.40 

17.34 

17.24 

83 

0700 

18.80 

17.63 

17.59 

17.55 

17.44 

83 

0720 

17.82 

17.76 

17.75 

17.62 

32 

0740 

19.40 

17.55 

17.71 

17.70 

17.53 

83 

1350 

20.4 

18.46 

18.36 

82 

1410 

21.8 

18.58 

18.77 

18.58 

79 

1620 

21.1 

18.7 

18.59 

19.12 

18.71 

79 

1650 

21.07 

18.69 

18.44 

18.83 

18.50 

79 

1710 

20.99 

18.56 

18.35 

18.58 

18.27 

79 

1730 

18.25 

18.43 

18.13 

79 
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Time 

■^1 

■^2 

^3 

^4 

H 

7/19 

1940 

18.5 

18.47 

18.42 

18.31 

86 

2000 

18.76 

18.03 

18.17 

84 

2040 

17.85 

17.79 

17.70 

17.54 

87 

2120 

17.76 

17.71 

17.64 

17.54 

87 

2140 

19.88 

17.77 

17.70 

17.65 

17.55 

87 

2200 

17.78 

17.60 

87 

7/20 

0700 

18.70 

17.41 

17.32 

17.29 

17.03 

86 

0740 

19.2 

17.51 

17.47 

17.27 

86 

0840 

19.22 

17.76 

17.70 

17.72 

17.48 

85 

0900 

19.3 

18.02 

18.02 

17.98 

17.84 

85 

0920 

19.3 

18.08 

18.07 

18.19 

17.88 

85 

1020 

19.63 

18.33 

18.33 

18.40 

18.12 

84 

1040 

19.92 

18.77 

18.92 

18.61 

80 

1240 

20.2 

19.1 

19.19 

19.25 

19.0 

78 

1300 

19.81 

19.12 

19.24 

19.22 

19.03 

79 

1320 

19.67 

19.08 

19.22 

19.24 

19.04 

88 

1800 

18.24 

18.80 

18.84 

19.16 

18.80 

84 

1900 

17.80 

18.86 

18.78 

18.87 

18.28 

83 

1920 

17.70 

18.77 

18.70 

18.79 

18.35 

84 

1940 

18.40 

18.70 

18.65 

18.28 

84 

2000 

18.30 

18.32 

18.38 

18.21 

85 

2020 

18.40 

13.27 

18.18 

17.74 

86 

2040 

18.14 

18.06 

17.85 

87 

2120 

18.24 

18.11 

18.08 

18.00 

17.83 

88 

2140 

17.97 

17.90 

17.71 

89 

2200 

18.63 

17.91 

17.84 

17.79 

17.64 

89 

2220 

19.00 

17.89 

17.83 

17.62 

90 

2230 

17.71 

17.71 

17.55 

91 

2300 

18.23 

17.38 

17.39 

17.23 

91 

2400 

17.20 

16.59 

16.59 

16.58 

16.57 

94 
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Time 

^1 

^2 

^3 

■^4 

H 

7/21 

0040 

16.9 

16.39 

16.37 

16.32 

16.20 

94 

0100 

16.6 

16.01 

16.04 

15.88 

93 

0120 

16.1 

15.83 

15.67 

15.66 

15.58 

96 

0405 

17.7 

16.42 

16.47 

16.38 

16.24 

98 

0425 

16.62 

16.61 

16.52 

16.40 

97 

0445 

18.4 

17.04 

17.08 

16.97 

16.83 

96 

0505 

18.4 

17.19 

17.20 

17.13 

17.07 

94 

0545 

18.2 

17.55 

17.57 

17.51 

17.39 

91 

0605 

17.67 

17.65 

17.59 

17.41 

89 

0625 

16.75 

16.69 

89 

0645 

18.3 

17.59 

17.56 

17.53 

17.34 

89 

0705 

18.19 

17.52 

17.48 

17.33 

89 

0725 

17.60 

17.56 

17.43 

89 

0745 

17.87 

17.79 

17.77 

17.63 

89 

0805 

17.73 

17.75 

17.64 

89 

0825 

17.71 

17.72 

17.57 

90 

0845 

19.04 

17.93 

17.93 

17.73 

91 

0905 

17.99 

18.03 

17.76 

89 

0945 

18.75 

17.75 

17.67 

17.75 

17.49 

89 

1005 

18.16 

17.63 

17.53 

17.73 

17.44 

88 

1025 

18.46 

17.60 

17.92 

17.58 

88 

1045 

18.4 

17.76 

17.63 

17.97 

17.60 

88 

1105 

17.72 

17.44 

17.42 

17.79 

17.39 

89 

1305 

17.66 

17.51 

17.94 

17.66 

90 

1325 

17.7 

17.72 

17.52 

17.85 

17.46 

90 

1345 

17.86 

17.97 

17.80 

17.97 

17.46 

90 

1405 

18.23 

18.06 

18.05 

18.14 

17.69 

90 

1505 

18.9 

18.53 

18.29 

18.60 

18.22 

88 

1620 

18.40 

18.79 

18.34 

86 

1720 

18.72 

18.51 

18.26 

18.55 

18.03 

85 

1945 

19.93 

19.93 

18.69 

18.86 

18.63 

79 

2030 

19.84 

18.45 

18.39 

18.31 

18.18 

85 

A-94 


Time 

^1 

^2 

■^3 

■^4 

H 

7/21 

2110 

18.50 

18.42 

18.32 

84 

2130 

19.5 

18.49 

18.46 

18.40 

18.29 

85 

2150 

18.59 

18.52 

18.41 

84 

2250 

18.28 

18.21 

18.15 

89 

2350 

17.95 

17.95 

17.89 

91 

7/22 


0130 

18.14 

18.11 

18.03 

92 

0150 

17.2 

18.02 

18.02 

18.01 

17.93 

91 

0450 

17.4 

17.41 

17.29 

17.31 

17.27 

92 

0550 

17.24 

17.21 

17.12 

93 

0610 

17.22 

16.83 

17.03 

16.98 

16.88 

94 

0710 

16.61 

16.60 

16.48 

96 

0730 

16.70 

16.70 

16.57 

97 

0750 

18.86 

16.67 

16.65 

16.68 

16.54 

97 

0810 

16.85 

16.86 

16.68 

96 

0830 

16.72 

16.64 

16.74 

16.58 

96 

0910 

16.66 

16.54 

16.65 

16.52 

97 

0930 

17.3 

16.74 

16.63 

16.77 

16.59 

97 

0950 

16.53 

16.72 

16.49 

97 

1010 

18.02 

16.75 

16.91 

16.66 

97 

1030 

18.51 

17.15 

17.18 

17.33 

17.05 

96 

1050 

18.6 

17.54 

17.64 

17.81 

17.51 

94 

1210 

17.24 

17.13 

17.27 

17.01 

90 

1230 

17.19 

17.41 

17.18 

91 

1333 

17.83 

17.80 

17.22 

92 

1350 

18.01 

17.92 

18.15 

17.59 

91 

1610 

19.81 

19.64 

19.90 

19.40 

87 

1630 

19.77 

20.51 

19.94 

87 

1650 

19.58 

20.35 

19.83 

86 
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Time 

^2 

»3 

*4 

H 

7/22 

1710 

19.68 

20.20 

19.85 

85 

1910 

19.34 

19.19 

19.29 

19.15 

90 

1950 

19.9 

18.78 

18.73 

18.53 

91 

2130 

17.97 

17.87 

17.82 

17.60 

93 

2150 

17.78 

17.64 

17.58 

17.33 

94 

7/23 

0310 

18.19 

17.83 

94 

0330 

18.17 

18.11 

17.85 

94 

0350 

18.24 

18.19 

17.91 

94 

0410 

18.23 

18.16 

17.99 

94 

0430 

18.25 

18.19 

18.15 

17.93 

94 

0450 

18.44 

18.39 

' 18.31 

18.09 

94 

0630 

18.56 

18.55 

18.49 

18.28 

91 

0730 

18.40 

18.60 

18.61 

18.58 

18.30 

90 

0750 

18.67 

18.65 

18.42 

89 

0810 

18.9 

18.80 

18.79 

18.58 

89 

0910 

19.06 

19.04 

18.97 

18.76 

88 

0930 

19.05 

18.96 

18.70 

88 

1410 

18.2 

19.25 

19.36 

18.95 

86 

1430 

18.2 

19.40 

19.33 

19.22 

86 

1440 

18.23 

19.51 

19.38 

19.58 

19.15 

87 

1450 

18.49 

19.69 

19.58 

19.75 

19.35 

86 

1505 

19.69 

19.83 

19.46 

85 

1515 

19.79 

19.91 

19.50 

84 

1525 

19.1 

20.02 

19.86 

20.10 

19.57 

64 

1535 

19.83 

19.97 

19.49 

63 

1645 

20.36 

19.87 

20.29 

19.93 

83 

1655 

19.83 

19.63 

19.90 

19.51 

83 

1705 

19.62 

19.47 

19.66 

19.34 

83 

1715 

19.33 

19.55 

19.21 

84 
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Time 

■^1 

^2 

^3 

•^4 

H 

7/23  ■ 

1725 

19.2 

19.23 

19.18 

19.43 

19.12 

85 

1735 

19.17 

19.18 

• 19.01 

19.28 

18.95 

85 

1745 

19.10 

18.93 

19.11 

18.82 

87 

1755 

19.03 

18.97 

18.80 

18.95 

18.65 

88 

1805 

19.2 

18.86 

18.74 

18.87 

18.59 

88 

2020 

17.7 

18.71 

18.68 

18.66 

18.31 

89 

2040 

17.7 

18.67 

18.64 

18.27 

89 

2340 

18.79 

18.78 

18.82 

18.53 

90 

2400 

18.4 

19.05 

18.91 

18.89 

18.60 

90 

7/24 


0020 

19.09 

19.06 

19.04 

18.86 

91 

0040 

19.12 

19.10 

18.95 

91 

0100 

19.19 

19.16 

19.15 

18.99 

90 

0120 

18.70 

19.08 

19.15 

19.15 

18.98 

90 

0240 

19.18 

19.13 

18.95 

87 

0300 

19.1 

19.17 

19.15 

18.98 

86 

0320 

19.04 

18.99 

19.00 

18.84 

87 

0420 

18.97 

18.93 

18.92 

18.76 

88 

0440 

18.81 

18.75 

18.75 

18.61 

88 

0500 

18.72 

18.70 

18.57 

88 

0520 

18.7 

18.81 

18.77 

18.73 

18.59 

89 

0540 

18.89 

18.83 

18.68 

89 

0600 

19.2 

18.98 

18.95 

18.91 

18.74 

88 

0620 

19.05 

19.00 

18.81 

88 

0640 

19.01 

18.99 

18.95 

18.77 

87 

0700 

18.94 

18.94 

18.74 

86 

0820 

19.03 

18.96 

19.26 

18.91 

81 

0840 

19.22 

19.07 

19.40 

19.07 

79 

0900 

19.15 

19.63 

19.23 

77 

0920 

19.14 

19.47 

19.18 

77 

0940 

19.32 

19.07 

19.46 

19.05 

78 
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Time 

^1 

■^2 

■^4 

H 

7/24 

1000 

19.39 

19.26 

19.66 

19.26 

78 

1020 

19.15 

19.28 

20.01 

19.41 

79 

1040 

19.32 

19.92 

19.35 

80 

1100 

18.75 

19.58 

19.87 

19.45 

80 

7/25 

2220 

17.7 

19.77 

19.82 

19.82 

19.27 

83 

2320 

17.86 

19.45 

19.46 

19.43 

19.13 

84 

7/26 

0020 

17.71 

19.39 

19.40 

19.37 

18.71 

85 

0120 

18.22 

19.05 

19.07 

19.13 

18.80 

87 

0320 

17.7 

18.96 

18.98 

18.98 

18.71 

89 

0420 

18.1 

18.86 

18.82 

18.84 

18.59 

90 

0440 

18.28 

18.41 

18.25 

91 

0500 

17.99 

18.26 

18.12 

89 

0600 

17.2 

18.23 

18.38 

18.42 

18.09 

91 

0620 

17.8 

18.17 

18.51 

18.21 

93 

0640 

19.05 

19.01 

18.60 

93 

0700 

19.76 

19.91 

18.58 

87 

0715 

19.92 

20.18 

19.76 

86 

0725 

19.79 

20.27 

19.76 

86 

0805 

19.39 

19.86 

19.42 

85 

0815 

19.33 

19.81 

19.38 

83 

0830 

19.43 

19.89 

19.48 

81 

0840 

19.32 

19.85 

19.38 

81 

0850 

19.26 

19.78 

19.29 

81 

0940 

18.86 

19.01 

18.68 

83 

1005 

18.83 

19.07 

18.74 

84 

1015 

18.77 

19.03 

18.70 

84 

1030 

18.80 

18.96 

18.71 

84 

1040 

18.98 

18.99 

18.81 

83 

1050 

19.11 

19.02 

18.96 

81 
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Time 


T 


4 


H 


7/26 


1105 

19.12 

19.06 

19.02 

81 

1115 

19.16 

19.04 

18.92 

80 

1130 

19.23 

19.10 

18.88 

78 

1140 

19.22 

19.10 

18.83 

77 

1155 

19.11 

19.05 

18.79 

78 

1210 

18.91 

18.88 

18.68 

80 

1235 

18.60 

18.58 

18.35 

80 

1245 

18.57 

18.52 

18.29 

80 

1300 

18.55 

18.53 

18.29 

78 

1310 

18.49' 

18.54 

18.20 

78 

1320 

18.47 

18.57 

18.22 

77 

1335 

18.37 

18.53 

18.20 

77 

1350 

18.12 

18.31 

17.95 

76 

1410 

17.56 

17.80 

17.40 

76 

1420 

77 

1435 

16.97 

17.18 

16.87 

78 

1450 

16.92 

17.16 

16.83 

79 

1500 

16.95 

17.19 

16.81 

78 

1515 

16.81 

17,09 

16.69 

78 

1525 

16.89 

17.27 

16.85 

77 

1540 

17.61 

17.17 

76 

1550 

17.32 

17.92 

17.43 

75 

APPENDIX  B:  CONPOTED  RESULTS 


This  section  Includes  cxaputed  values  of  Richardson's  number,  the 
momentum  and  heat  flux,  dlffusivity,  and  the  dissipation  rate  and  twqpera- 
ture  structure  function,  nie  mean  wind  spaed  is  Included  for  convenience. 
Figures  showing  the  ten{>oral  variation  of  Richardson's  number  and  momentum 
fliix,  and  the  dependence  of  nnmentum  flux  on  wind  speed  are  also  included. 
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Computed  Results 
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APPENDIX  C:  AEROSOL  RESULTS 


This  section  is  excerpted  from  a Masters  Thesis  written  by  LT  Alan 
Simoncek,  USN.  It  includes  data  and  analysis  for  the  cruise  reported 
here  and  also  for  data  obtained  off  the  Florida  Coast.  Both  are  included 
for  comparison  purposes.  All  of  the  data  was  gathered  by  personnel  from 
Calspan  Corporation,  Buffalo,  New  York. 
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I.  INTRODUCTION 


The  military  is  currently  very  interested  in  the  perform- 
ance of  electro-optical  weapons  systems  in  an  atmosphere  of 
varying  turbidity.  For  example,  a number  of  electro-optical 
systems  which  utilize  the  visible  as  well  as  IR  wavelengths 
are  being  developed  by  the  Navy  for  use  in  surveillance  and 
intelligence  gathering  operations  in  the  marine  boundary 
layer.  These  systems  are  limited  by  the  extinction  of  the 
propagated  energy  due  to  absorption  and  scattering  by  aerosols. 
The  effect  of  absorption  depends  on  the  composition  of  the 
particulates  and  wavelength  of  the  energy  and  the  effect  of 
scattering  depends  on  the  concentration  and  size  of  the  scat- 
terers.  For  most  applications  the  scattering  processes  in 
the  atmosphere  are  caused  by  particles  of  size  comparable  to 
the  wavelength  of  the  radiation. 

The  size  distribution  of  the  marine  aerosol  is  known  to 
depend  upon  the  wind  speed,  relative  humidity,  stability,  and 
air  mass  trajectory.  In  order  to  evaluate  accurately  and  pre- 
dict the  atmospheric  effects  on  these  electro-optic  systems, 
it  is  necessary  to  know  the  dependence  of  the  aerosol  size 
distribution  on  the  foregoing  meteorological  parameters. 

The  nature  of  the  aerosol  size  distribution  in  a coastal 
marine  environment  is  investigated  in  this  study.  Data  from 
aerosol  observations  off  the  coast  of  Panama  City,  Florida 
and  off  the  Southern  California  coast  near  the  Channel  Islands 
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were  analyzed.  These  coastal  regimes,  which  represent  a 
mixture  of  continental  and  marine  aerosols , should  contain 
aerosol  distributions  somewhat  different  from  the  typical 
marine  environment.  The  relationship  of  the  coastal  marine 
aerosol  to  wind  speed,  relative  humidity,  stability,  and 
sub-synoptic  circulation  is  examined.  Furthermore,  an 
attempt  is  made  to  evaluate  the  use  of  the  friction  veloc- 
ity as  a valid  aerosol  distribution  predictor. 
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II.  BACKGROUND 


A.  THE  ATMOSPHERIC  AEROSOL 

With  the  recent  increasing  concern  over  the  pollution 
of  our  atmospheric  environment,  the  examination  of  the 
tropospheric  aerosols  has  also  increased.  Particulate  mat- 
ter enters  the  atmosphere  through  either  natural  or  man-made 
processes;  approximately  10%  of  the  total  concentration  is 
believed  to  originate  from  combustion  and  industrial  pro- 
cesses while  the  natural  sources,  including  soil  dust, 
volcanoes,  and  oceans  account  for  the  remaining  90%.  The 
size  range  of  aerosols  observed  by  current  methods  extends 
from  lO”  y to  10  y radius  (1  y = 10  m = micron).  Depend- 
ing on  their  size,  amount  of  soluble  matter,  and  the  rela- 
tive humidity,  these  particles  may  act  as  condensation 
nuclei  and  aid  in  the  precipitation  process. 

Mason  (1975)  classified  condensation  nuclei  into  three 
groups  according  to  radius:  Aitken  (<  0.1  y).  Large  (0.1  y 
- 1 y),  and  Giant  (>  1 y)  particles.  Essentially,  Aitken 
nuclei  are  produced  by  man-made  sources  and  larger  nuclei 
by  natural  processes.  Therefore,  it  is  not  surprising  to 
see  Aitken  nuclei  dominate  the  size  distribution  spectrum 
over  continents.  The  marine  aerosol  above  .1  y is  composed 
of  sea-salt  particles  produced  by  spray  and  bubble  bursting 
mechanisms  on  the  water  surface.  These  mechanisms  are  quite 
complex  and  their  contribution  to  the  size  distribution  will 
be  discussed  in  detail  later. 
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The  vertical  profiles  of  trace  constituents  that  are  pro- 
duced over  the  continents  have  been  shown  to  be  rather  uni- 
form in  space  and  time  above  5 km.  This  "background"  aerosol 
is  affected  slightly  by  anthropogenic  activities  and  is  far 
away  from  local  natural  sources.  Its  number  concentration 
is  almost  identical  with  the  concentration  of  Aitken  nuclei 

over  ocean  areas.  Past  experiments  have  shown  the  concentra- 

- 3 

tion  of  the  background  aerosol  to  be  about  300  cm  over 
remote  ocean  areas.  However,  a recent  experimental  cruise 

(R/V  Meteor)  by  Junge  and  Jaenicke,  1971  in  the  mid  Atlantic 

. - 3 

yielded  observed  concentrations  of  600  cm  . Measurements 

by  Hidy,  ejt  al.  (1973)  on  San  Nicolas  Island,  130  km,  west- 

southwest  of  Los  Angeles,  have  shown  the  background  aerosol 

to  be  a mixture  of  material  from  both  marine  and  continental 

. - 3 

sources  with  an  average  concentration  of  2400  cm  . Samples 
taken  over  oceans  of  thd  South  Atlantic  (Meszaros  and  Vissy, 
1974)  resulted  in  Aitken  particle  counts  of  between  300- 
450  cm” ^ . 

B.  CHARACTERISTICS  OF  THE  MARINE  AEROSOL 

An  idealized  size  distribution  of  the  continental  and 
marine  aerosol  is  supplied  by  Junge  (1972)  in  Figure  1. 

The  significant  feature  is  the  shift  of  the  maximum  of  par- 
ticles as  a function  of  total  particle  concentration.  Over 
the  ocean  the  sea-salt  aerosol,  which  is  usually  confined 
to  the  lower  2 km,  is  superimposed  on  the  background  aerosol. 
Junge  reasoned  that  the  concentration  of  the  background 
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Figure  1,  Idealized  Size  Diatributlona  of ‘Continental' 
and  Marine  Aerosols 
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cloud  nuclei  decreases  below  the  marine  inversion  due  to 


the  effect  of  washout,  or  coalescence,  due  to  the  larger 
water  droplets  in  the  cumulus  clouds. 

Another  aspect  of  the  aerosol  distribution  is  the 
slope  of  the  number  density  versus  radius  curve.  Fried- 
lander  (1961)  proposed  a theory  of  self-preserving  size 
distributions  which  helps  to  explain  why  all  atmospheric 
size  distributions  are  similar.  He  proposed  that  the 
similarities  can  be  explained  by  solutions  to  the  kinetic 
equation  which  describe  the  relationship  between  particle 
size  distribution  and  time.  Experimental  results  have  in- 
dicated that  the  size  distribution  over  a particular  range 
of  sizes  of  continental  aerosol  has  a -4  slope  and  follows 
the  relation 


dN 

H?  = 


C<))  r 


(1) 


3 

where  N is  the  number  of  particles/cm  , r the  particle 
radius,  C a constant,  and  (|)  the  volume  of  particles  per 
unit  volume  of  aerosol. 

Blifford  (1970)  measured  the  size  and  number  distribu- 
tions of  atmospheric  aerosols  at  various  altitudes  over  the 
ocean  250  km  west  of  Santa  Barbara,  California.  Samples 
were  taken  by  an  aircraft  equipped  with  a jet  impactor  and 
the  data  was  obtained  from  direct  microscopic  counting  tech- 
niques in  the  laboratory.  The  aerosol  distribution  at  ap- 
proximately 15  meters  above  the  sea  surface  is  presented  in 
Figure  2.  The  curve  has  a rather  steep  negative  slope  at 
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RADIUS  (MCRONS) 


Figure  2,  Size  Distribution  at  15  Meters,  250  km  West 
of  Santa  Barbara 
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the  small  particle  end  which  becomes  slightly  positive  at 
around  .4  y radius.  For  particles  larger  than  .8  y,  a 
fairly  constant  slope  of  about  -2  to  -3  is  observed. 

The  results  of  the  R/V  Meteor  experiment,  where  several 
aerosol  counters  were  used,  are  shown  in  Figure  3.  Above 
10  y the  exponent  of  a power  function  fit  to  the  data  is 
approximately  -6  and  between  0.3  y and  10  y it  is  variable 
but  on  the  average  around  -3.  The  maximum  of  the  size  dis- 
tribution occurred  at  0.3  y with  a secondary  maximum  at 
0.03  y . 

It  is  possible  that,  due  to  increased  human  activity 
in  the  Northern  Hemisphere,  Junge  and  Jaenicke's  Atlantic 
experiment  did  not  explore  the  undisturbed  marine  environ- 
ment. Meszaros  and  Vissy  (1974)  describe  the  results  of 
aerosol  samples  taken  over  the  oceans  of  the  Southern  Hemi- 
sphere by  means  of  membrane  filters.  An  example  of  the 
number  concentration  and  size  distribution  over  the  Atlantic 
between  (a)  0°  and  20°  South  and  (b)  40°  and  60°  South  can 
be  found  in  Figure  4.  Chemical  analyses  were  performed 
and  it  was  observed  that  the  maxima  in  the  concentrations 
of  all  particles  and  of  sodium  chloride  particles  occur  at 
approximately  .1  y radius  in  both  cases.  Up  to  . 5 y radius 
the  slope  of  the  distribution  is  approximately  -5.  Between 
0.5  y and  1.5  y,  however,  the  decrease  of  the  concentration 
with  increasing  particle  size  is  very  moderate  (-1  to  -2), 
while  for  radii  larger  than  1,5  y the  slope  is  close  to  -3. 
This  has  been  interpreted  to  indicate  that  the  form  of  the 
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Figure  •3.  Results  of  the  R/V  Meter  Cruise.  Measurement 
Systems  Used;  o— o Combination  of  CCN  Counter, 

Optical  Counter,  and  Impactors; Double 

Stage  Impactor  (Junge  and  Jaenicke,  1971) 
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distribution  is  produced  by  the  combined  effect  of  particles 
formed  in  different  ways. 

Oceanic  measurements  have  shown  that  the  concentration 
of  sea-salt  particles  decreases  exponentially  with  height, 
with  little  variation  of  the  size  distribution.  Ericksson 
(1959)  reported  that  there  exists  a level  a few  hundred 
feet  above  the  surface  where  the  concentration  decreases 
with  height  in  periods  of  high  wind  force  and  increases  with 
height  in  lower  wind  forces.  He  reasoned  that  there  is 
little  or  no  production  of  sea-salt  in  regions  of  light 
winds  and  that  coagulation  and  fallout  in  the  lowe.  'cvels 
combined  with  horizontal  transport  due  to  vertical  shear 
produce  a maximum  concentration  at  some  upper  level. 

Toba  (1965a  and  b)  proposed  that  the  average  decrease 
in  concentration  with  height  can  be  explained  by  a combina- 
tion of  sedimentation,  diffusion,  convective  processes  and 
the  humidity  distribution.  He  suggested  that  the  line  be- 
tween the  aerosol  vertical  distribution  and  the  process  of 
production  of  sea-salt  particles  at  the  sea  surface  is  found 
within  the  lowest  layer  of  the  atmosphere  where  the  eddy 
diffusivity  and  relative  humidity  sharply  change. 

The  distribution  of  eddy  diffusivity  near  the  sea  surface 
is  closely  related  to  the  wind  speed.  The  larger  the  eddy 
diffusivity  near  the  surface  the  more  sea-salt  particles 
that  will  be  supplied.  Toba  considered  eddy  diffusivity  in 
the  form 


D = kU.(Z  + Z) 
» o 


(2) 
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where  k is  the  von  Karman  constant,  Z the  height  above 
the  sea  surface,  the  roughness  length,  and  the 

friction  velocity  which  is  a function  of  the  momentum  trans- 
fer over  the  sea. 

The  relative  humidity  in  the  first  few  meters  over  the 
ocean  is  known  to  decrease  rapidly  with  height.  The  par- 
ticles produced  near  the  surface  in  a region  of  high  humid- 
ity grow  larger  and  thus  have  a greater  terminal  velocity 
due  to  gravity  than  those  at  the  top  of  this  layer. 

During  light  winds  the  number  concentration  near  the  sea 
surface  increases  with  height.  Since  it  results  from  a non- 
steady state,  an  inversion  of  vertical  gradient  of  the  par- 
ticle concentration  is  most  likely  to  be  found  in  small  par- 
ticles which  have  a longer  residence  time.  Ericksson  (1959) 
computed  the  fall  velocities  for  given  relative  humidities, 
salinity,  and  radius.  During  high  wind  periods,  giant  size 
sea-salt  particles  are  produced  at  the  surface  and  through 
the  diffusion  process  are  mixed  throughout  the  atmosphere. 

The  largest  particles  may  fall  back  into  the  ocean  due  to 
excessive  terminal  velocity  or  be  entrained  in  the  wave 
crests.  Smaller  particles  are  free  to  rise  to  cloud  height 
where  coalescence  with  larger  cloud  drops  and  washout  usually 
occur. 

-14 

Measurements  of  salt  nuclei  greater  than  10  gm  over 
the  North  Atlantic  by  Moore  and  Mason  (1954)  revealed  the 
existence  of  two  distinct  types  of  size  distributions  (Type 
I and  II).  The  curves  for  the  observed  Type  I and  Type  II 
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nuclei  distributions  are  reproduced  in  Figure  5.  Type  I 
distributions  were  observed  for  wind  speeds  between  6-15 
m/sec  and  were  thought  to  be  residuals  of  spray  droplets 
produced  by  breaking  waves.  The  presence  of  a discontinuity 
or  a sharp  change  of  slope  in  the  Type  I distribution  was 
explained  in  terms  of  a loss  of  the  larger  nuclei  by  sedi- 
mentation. In  strong  winds,  the  part  of  the  curve  to  the 
right  of  the  discontinuity  probably  represents  a state  of 
equilibrium  between  production  and  loss  by  sedimentation. 

In  light  winds  and  stable  conditions  the  slope  should  be 
steeper  due  to  the  fact  that  the  loss  by  sedimentation  is 
greater  than  production  and  larger  nuclei  are  not  easily 
transported  vertically  under  stable  conditions.  The  Type  II 
distributions  were  only  observed  when  the  wind  speeds  were 
less  than  7 m/sec  and  resembled  a high  concentration  conti- 
nental aerosol.  In  winds  of  up  to  15  m/sec  the  measured 

concentrations  of  large  sea-salt  nuclei  rarely  exceeded  10 

-3 

cm 

The  effect  of  stability  on  the  concentration  of  atmos- 
pheric condensation  nuclei  was  well  documented  by  Moore 
(1952).  He  used  an  Aitken  counter  to  measure  the  relation- 
ship between  concentration  of  nuclei  and  the  intensity  of 
vertical  mixing  over  the  North  Atlantic.  The  results  indi- 
cated a decrease  by  as  much  as  a factor  of  4 in  the  number 
of  Aitken  particles  near  the  surface  on  days  with  cumulus 
clouds  as  compared  to  days  with  stratus  clouds.  This  would 
indicate  that  convection  plays  an  important  role,  at  least 
in  the  transport  of  smaller  particles. 
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Figure  5,  The  Size  Distribution  of  Moore  and  Mason’s 
Type  I and  Type  II  Nuclei 
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Chemical  analyses  by  various  investigators  have  indi- 
cated that  between  0.1  y and  1.0  y radius  the  marine 
aerosol  is  composed  of  a background  component  of  continental 
origin  and  a sea-salt  component.  Sodium  chloride  was  found 
to  predominate  above  1 y radius  while  particles  of  continen- 
tal origin  predominate  below  0.1  y radius.  Results  of  a 
cruise  off  the  Grand  Banks  in  the  North  Atlantic  (Ruskin, 
et  al. , 1976)  indicated  that  the  continental  particles  are 
composed  of  sulfate  compounds  and  a smaller  amount  of  sul- 
furic acid.  Aerosols  over  remote  ocean  areas  (Meszaros  and 
Vissy,  1974)  were  shown  to  be  comprised  of  variable  concen- 
trations of  ammonium  sulfate,  sulfuric  acid,  sodium  chloride, 
and  particles  similar  in  structure  to  ammonium  sulfate.  The 
sum  of  these  four  types  of  identified  particles  accounted 
for  75-95  percent  of  all  particles  greater  than  .3  y radius. 
In  other  words,  practically  all  the  particles  in  a pure  mar- 
ine atmosphere,  undisturbed  by  continental  particle  sources, 
are  soluble  in  water. 

C.  RELATIVE  HUMIDITY  EFFECTS 

A solid  particle  which  is  composed  wholly,  or  in  part, 
of  a pure  water-soluble  substance  will  undergo  a sudden 
transition  to  a saturated  solution  droplet  when  some  critical 
value  of  relative  humidity,  less  than  100%,  is  reached.  The 
relative  humidity  at  which  this  transition  occurs  depends  on 
the  size  and  chemical  composition  of  the  particle.  The 
smaller  the  particle,  the  lower  the  critical  humidity.  Below 
the  transition  point,  solid  particles  acquire  small  amounts 
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of  water’  by  the  process  of  adsorption.  At  relative  humidi- 
ties above  the  transition  point,  a particle  (or,  more  pro- 
perly, an  aqueous  solution  droplet)  grows  by  the  absorption 
of  water  vapor  (Fitzgerald,  1975). 

A pure  water  droplet  is  said  to  be  in  equilibrium  with 
its  surroundings  if  it  neither  evaporates  nor  grows.  This 
only  occurs  when  the  equilibrium  vapor  pressure  over  the 
surface  of  the  droplet  is  equal  to  the  vapor  pressure  of  the 
surrounding  air.  Winkler  (1973)  describes  the  equilibrium 
growth  of  aerosol  particles  due  to  humidity  as  complex  and 
depending  on  the  relative  proportion  of  soluble  and  insoluble 
material  in  the  particles  and  on  the  chemical  composition  of 
the  soluble  component.  Complex  ionic  mixtures,  similar  to 
those  present  in  a:  mospheric  aerosols,  show  material  in- 
fluences and  lower  the  water  vapor  pressure  to  a much  less 
degree  than  the  same  amount  of  pure  salts.  In  such  complex 
mixtures  the  various  salts  become  dissolved  only  gradually 
with  increasing  relative  humidity  until  at  a sufficiently 
high  humidity  all  soluble  material  is  in  solution. 

Measurements  have  shown  that  with  increasing  humidity  a 
sodium  chloride  crystal  undergoes  a phase  transition  to  a 
saturated  solution  droplet  at  a relative  humidity  of  approxi- 
mately 78%.  Figure  6 describes  how  the  equilibrium  radii  of 
droplets  containing  specified  masses  of  sodium  chloride  vary 
with  the  relative  humidity.  The  equilibrium  radius  of  the 
droplet  increases  with  increasing  humidity  until  the  air 
becomes  supersaturated  by  a critical  amount,  corresponding  to 
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Figure 


Equilibrium  Relative  Humidity  eind  Corresponding 
Radii  (Mason,  1975) 
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the  maximum  of  the  curve  in  this  figure.  If  this  supersatu- 
ration were  maintained,  theoretically  the  droplet  will  grow 
without  bound.  With  decreasing  humidity  a sodium  chloride 
solution  droplet  crystallizes  at  a humidity  between  35-45%. 
Since  the  relative  humidity  at  a height  of  about  15  meters 
over  the  ocean  surface  goes  below  40%  very  infrequently,  sea- 
salt  droplets  will  have  little  opportunity  to  crystallize 
(Fitzgerald  and  Ruskin,  1977). 

Since  the  later  discussion  refers  to  the  distribution 
of  sea-salt  particles  by  bolt  mass  weight  of  salt  (grams), 
radius  of  dry  crystals  (p)  or  radius  at  ambient  humidity 
(p),  the  scale  in  figure  7 is  furnished  as  a reference. 

D.  THE  PRODUCTION  OF  AIRHOLE  SEA-SALT 

Although  the  spectrum  of  the  marine  aerosol  above  .1  p 
radius  is  known  to  consist  of  sea-salt  particles,  very  little 
is  certain  about  the  concentration  and  mechanisms  of  produc- 
tion. Because  of  the  smallness  of  the  particles  and  limita- 
tions of  the  sampling  equipment,  earlier  experiments  did  not 
measure  the  quantity  of  sea-salt  particles  much  less  than 
10  gm. 

Woodcock  (1953)  determined  that  the  mass  distribution  of 
"giant"  (>  10  gm)  sea-salt  nuclei  varies  with  wind  speed. 
Increases  in  the  amount  of  air-borne  salt  near  cloud  bases 
were  shown  to  be  related  to  increases  in  wind  speed  at  the 
sea  surface,  with  the  greatest  proportionate  increase  in  par- 
ticle number  occurring  at  the  large  end  of  the  weight  range. 
The  results  of  Woodcock's  measurements  for  wind  forces  of  1, 
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3,  5,  and  7 on  the  Beaufort  scale  are  shown  in  Figure  7.  The 
line  (a)  gives  the  size  distribution  of  continental  aerosol 
for  comparison.  The  line  (b)  is  an  extrapolated  size  distri- 
bution of  the  marine  aerosol.  Chemical  analysis  of  Woodcock's 
bulk  aerosol  samples  between  .1  y and  1 y indicated  a maximum 
of  sea-salt  around  0.3  y and  a lower  limit  in  the  vicinity  of 
.1  y radius.  These  distributions  indicate  total  concentra- 
tions of  all  sea-salt  particles  of  no  higher  than  a few  per 
cubic  centimeter  (Junge,  1972).  According  to  Mason  (1975), 

over  a rough  sea  the  concentration  of  sea-salt  particles 

- 3 

greater  than  2 y radius  rarely  exceeds  1 cm  and  the  total 

. - 3 

concentration  of  all  salt  particles  rarely  exceeds  10  cm 

Moore  (1952)  observed  a distinct  correlation  between  wind 
speed  and  concentration  of  sea-salt  larger  than  lO”^^  gm  up 

to  wind  speeds  of  15  m/sec.  He  also  found  a linear  increase 

. . -9 

in  concentration  of  particles  larger  than  10  gm  with  in- 
crease in  v/ave  height.  Results  of  experiments  by  Monahan 
(1968)  reveal  an  abrupt  increase  in  concentration  of  sea 
water  droplets  larger  than  45  y radius  at  a wind  speed  of 
approximately  9 m/sec,  measured  47  cm  above  the  sea  surface. 

Moore  (1952)  also  analyzed  the  visibility  observations 
at  two  ocean  weather  ships  and  determined  that  the  opacity 
for  a given  humidity  increases  with  wind  speed.  He  attri- 
buted this  increase  to  an  observed  increase  in  the  concentra- 
tion of  large  nuclei.  Another  result  indicated  that  at  lower 
humidities,  the  increase  in  opacity  was  more  pronounced,  and 
Moore  believed  this  was  due  to  the  dehydration  of  larger 
droplets.  These  conclusions  would  indicate  that  the  aerosol 
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Figure  7.  Size  Scale  and  Average  Size  Distribution 

of  Sea-Salt  Nuclei  Measured  by  A.H.  Woodcock 
(Mason,  1975) 


distribution  is  more  variable  and  sensitive  to  wind  speed 
in  drier  air,  and  this  feature  should  be  most  noticeable  in 
the  larger  size  ranges. 

As  the  wind  speed  increases  over  the  ocean,  gravity 
waves  are  generated  and  begin  to  break  at  a critical  wind 
speed  generally  agreed  upon  to  be  near  7 m/sec.  Air  that 
is  entrained  by  these  breaking  "whitecaps"  rises  to  the  sur- 
face sometime  later  in  the  form  of  bubbles.  The  principal 
mechanisms  of  sea-salt  production  are  thought  to  be  the 
direct  spraying  of  droplets  off  the  crests  of  breaking  waves 
and  the  bursting  of  bubbles  in  areas  of  whitecaps  and  foam. 
Droplets  produced  by  direct  spraying  are  generally  larger 
than  45  y and,  due  to  large  fall  velocities,  are  not  air- 
borne long  enough  to  evaporate  and  become  light  enough  to 
be  transported  upward  (Monahan,  1968).  Toba's  model  (1965b) 
showed  that  the  net  production  of  sea-salt  particles  at  the 
sea  surface  seems  to  increase  with  particle  mass  even  beyond 
lO”  gm  (20  y),  but  that  the  transport  by  eddy  diffusion  is 

not  sufficient  to  carry  the  particles  upward  against  gravity 

— 8 

beyond  this  size.  The  presence  of  particles  larger  than  10 
gm  in  the  atmosphere  is  generally  attributed  to  coalescence 
of  sea-salt  droplets  within  and  below  clouds. 

Some  examples  of  residence  times  for  different  sea-salt 
particle  sizes  taken  from  Junge  (1972)  are  found  in  Table  I. 
It  would  seem  then  that  particles  in  the  .1  y - 20  y range, 
at  least,  are  produced  by  the  bursting  of  bubbles. 

In  efforts  to  photograph  the  rupture  of  the  surface 
bubble  film,  Kientzler,  et  al.  (1954)  found  that  their  camera 
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exposure  was  too  long  to  capture  this  rapid  phenomenon.  How- 
ever, they  V7ere  able  to  see  the  formation  of  the  "Rayleigh" 
jet  which  projects  upward,  continues  to  rise  as  a thin  column, 
and  then  breaks  into  droplets  of  varying  sizes.  Day  (1964) 
describes  this  process  in  the  following  manner.  Each  bubble, 
as  it  reaches  the  surface,  develops  a spherical  film-cap 
which  drains,  thins,  and  bursts.  Fragments  of  the  film  are 
thrown  out  and  are  dragged  upward  by  the  air  which  escapes 
from  the  bubble  orifice.  Water,  rushing  down  the  sides  of 
the  bubble  cavity,  emerges  from  the  center  as  a narrow  jet. 

A schematic  of  this  process  is  shown  in  Figure  8.  The  larger 
drops  (L)  are  formed  by  disintegration  of  the  jet  (J). 

Smaller  particles  (S)  are  formed  by  bursting  of  the  bubble 
film. 

Kientzler's  experiment  was  significant  in  that  no  drop- 
lets of  large  enough  size  to  be  resolved  by  the  film  and 
optical  system  were  observed  from  .2  - 1.8  mm  diameter  bubbles 
until  after  the  jet  formation.  This  was  interpreted  to  indi- 
cate that  the  larger  droplets  are  not  produced  when  the 
bubble  film  is  broken.  On  the  average,  the  droplets  produced 
by  the  jet  mechanism  were  approximately  1/10  of  the  original 
bubble  size.  1 mm  diameter  bubbles  were  observed  to  produce 
droplets  of  approximately  50  p radius.  The  smallest  observed 
were  of  2 p radius  and  deduced  to  have  been  formed  by  a 
bubble  of  approximately  .04  mm  diameter.  Therefore,  the  jet 
mechanism  can  be  considered  a source  of  salt  particles  greater 
than  10”^^  gm  (1  p). 
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Figure  8,  The  Formation  of  Sea-Salt  Droplets  by  the 
Bursting  of  Bubbles  (Mason,  1975) 
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Mason  (1954)  utilized  a cloud  chamber  to  study  bubble 
behavior  in  both  distilled  and  salt  water.  After  expansion, 
a dense  cloud  of  tiny  droplets  was  observed  rising  vertical- 
ly in  the  space  above  the  salt  water,  but  not  above  the  dis- 
tilled water.  Bubbles  of  3 mm  diameter  produced  100-200  of 

these  condensation  nuclei,  the  majority  of  which  are  esti- 

-15  -14 

mated  to  have  salt  contents  between  10  gm  and  2 x 10  gm. 
This  v;ould  correspond  to  droplets  of  approximately  .1  y-  to 
.3  y radius  at  80%  relative  humidity.  Mason  also  observed 
a second  group  of  droplets  produced  by  the  shattering  of  the 
bubble  film.  These  were  projected  sideways  at  an  angle  of 
ten  to  15  degrees  above  the  horizontal  and  slightly  larger, 
containing  between  2 x 10~^^  - 5 x 10~^*^  gms  of  salt.  How- 
ever, the  numbers  of  these  droplets  were  always  small,  on  the 
average,  there  was  only  about  one  droplet  in  this  size  range. 

The  number  of  droplets  which  rise  vertically  from  a 
bursting  bubble  is  strongly  I’elated  to  the  state  of  compres- 
sion of  the  film  of  organic  material  on  the  water  surface. 
Paterson  and  Spillane  (1969)  have  shown  that  with  an  increase 
of  film  pressure  the  number  of  nuclei  produced  decreases 
markedly.  This  would  indicate  that  the  production  of  sea- 
salt  droplets  originating  from  the  bubble  film  mechanism  would 
be  suppressed  in  regions  of  high  organic  activity  on  the  sea 
surface.  Aerosol  samples  taken  by  Woodcock  (1972)  over 
Hawaiian  and  Alaskan  seas  may  help  explain  where  the  transi- 
tion between  the  jet  and  film  sea-salt  production  mechanisms 
occurs.  His  observations,  using  an  improved  slide  collection 


C-32 


technique,  show  an  increased  average  particle  proauction  for 

- 14 

sea  salt  particles  less  than  2 x 10  gra  ( . 3 p radius)  in 
Hawaii  where  marine  organic  productivity  is  low.  In  contrast, 
the  mean  distribution  curve  for  particles  over  the  organical- 
ly rich  Gulf  of  Alaska  fails  to  indicate  an  increased  slope 
of  the  concentration  curve  among  particles  of  the  same  size 
range.  These  curves  are  shown  in  Figure  9.  The  presence  of 
surface  active  films  arising  from  the  biologically  productive 
Alaskan  waters  is  thought  to  suppress  the  production  of  film 
droplets . 

Statistical  analysis  by  Meszaros  and  Vissy  (1974)  showed 
that  with  increasing  particle  radius  the  correlation  between 
wind  speed  and  chloride  concentration  increased.  This  meant 
that  smaller  chloride  particles  are  formed  by  the  bubble 
film  mechanism  than  by  direct  spraying.  The  distribution 
curve  gives  evidence  that  the  transition  between  these  two 
chloride  formation  mechanisms  lies  between  .2  p and  .4  p. 

Thus  the  maximum  at  . 1 p gives  the  maximum  of  chloride  par- 
ticles formed  by  the  bubble  film  process. 

Moore  (1952)  found  evidence  that  the  particle  concentra- 
tions below  1 p are  not  correlated  with  wind  speed.  This 
would  indicate  that  most  of  the  particles  between  .1  p and 
1 p are  not  produced  by  bursting  bubbles.  Other  experiments 
using  the  effects  of  relative  humidity  on  particle  growth 
indicate  that  a considerable  proportion  of  marine  particles 
between  .1  p and  1 p must  differ  in  composition  from  sea- 
salt  (Junge,  1972).  Meszaros  and  Vissy  (1974)  found  that, 
in  this  size  range,  sodium  chloride  varied  from  4-50%  of  the 
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Figure  9,  Distribution  Curves  over  Alaskan  and  Hawaiian 
Waters  (Woodcock,  1972) 


9 


C-3 


concentration  for  all  particles.  The  observations  by  Kidy 
et  (1974)  off  the  coast  of  Southern  California  revealed 

that  11%  of  the  aerosol  sampled  contained  sea-salt,  the 
remainder  being  a combination  of  sulfates,  nitrates  and  soil 
dust . 

E.  AEROSOL  MODEL 

Recently,  various  aerosol  models  have  been  developed  in 
an  attempt  to  accurately  describe  marine  aerosol  distribu- 
tions as  a function  of  one  or  more  parameters.  This  is  es- 
sential for  the  calculation  of  optical  propagation  through 
the  atmosphere  as  aerosols  scatter  and  absorb  energy.  Since 
the  aerosol  distribution  is  known  to  be  dependent  on  relative 
humidity  and  wind  speed,  these  two  variables  usually  are  the 
key  parameters  of  each  model. 

One  model  in  particular  has  been  developed  by  Fitzgerald 
and  Ruskin  (1977)  on  the  basis  of  the  North  Atlantic  observa- 
tions. They  applied  the  effects  of  relative  humidity  on  the 
equilibrium  growth  of  aerosol  particles  to  the  sea-salt  mass 
distribution  determined  by  Lovett  (1975)  in  the  North 
Atlantic.  Lovett  presents  empirical  log  radius  mass  distri- 
butions in  the  form  of  the  following  power  law: 

^ = C r/'"  (3) 

d log  r^ 

where  r^  is  the  dry  particle  radius  and  C and  v depend 
on  the  wind  speed  V in  the  following  manner: 

V = 3.317  - .03  V (4) 
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and  C = 0.2  - 0.0196  V + 0.0121  (5) 

These  expressions  are  valid  only  over  a wind  speed  range  of 
3-17  m/sec~^. 

Formulae  have  been  derived  (Fitzgerald,  1975)  for  the 
equilibrium  size  of  aerosol  particles  composed  of  a single 
pure  salt  as  a function  of  relative  humidity.  For  a sodium 
chloride  particle  the  relationship  between  particle  radius 
and  relative  humidity  may  be  expressed  as 

r = a r^®  (6) 

where  a and  6 are  functions  of  the  relative  humidity  as 
described  by  Fitzgerald  (1975).  Equations  (3)  and  (6)  are 
combined  to  describe  the  aerosol  size  distribution  as  a 
function  of  relative  humidity  and  wind  speed,  giving 


— SN = £ („'''«)  . 

d log  r 6 


(7) 


Comparison  between  the  aerosol  distributions  derived  from 
the  above  model  and  those  observed  in  two  coastal  marine 
environments  is  made  within  this  study. 
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III.  TURBULENCE  THEORY 

A.  BOUNDARY  LAYER  CONSIDERATIONS 

I The  importance  of  turbulent  exchange  processes  in  the 

surface  boundary  layer  has  long  been  recognized.  Panofsky 
(1969)  describes  atmospheric  turbulence  as  consisting  of 
horizon!  il  and  vertical  eddies  by  which  the  air  is  mixed. 

The  two  mechanisms  by  which  eddies  are  formed  in  the  atmos- 
phere are  heating  from  below  and  wind  shear.  Heating  pro- 
duces convection  and  the  change  in  wind  speed  with  height 
produces  mechanical  turbulence.  Because  there  is  no  wind 
at  ground  level,  and  there  is  usually  some  wind  above  the 
ground,  mechanical  turbulence  is  common.  This  type  of  tur- 
bulence increases  with  increasing  wind  speed  (at  a given 
height)  and  is  greater  over  rough  terrain  than  over  smooth 
terrain.  The  terrain  roughness  is  usually  characterized  by 
a roughness  length,  , which  is  proportional  to  the  size 

of  the  eddies  that  can  exist.  The  relative  importance  of 
heat  convection  and  mechanical  turbulence  is  characterized 
by  the  Richardson  number,  R^  . The  Richardson  number  is  a 
measure  of  the  relative  rate  of  conversion  of  convective  to 
mechanical  energy.  For  example,  negative  Richardson  numbers 
of  large  magnitude  indicate  that  convection  predominates  re- 
sulting in  strong  vertical  motion.  As  the  mechanical  turbu- 
lence increases,  the  Richardson  number  approaches  zero. 
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Finally,  as  the  Richardson  number  becomes  positive,  the  ther- 
mal stratification  becomes  stable  and  damps  the  mechanical 
turbulence.  For  > 0.25,  vertical  mixing  disappears. 

The  effect  of  the  wind  on  the  underlying  surface  is 
termed  the  shearing  or  Reynolds  stress , t , and  is  character- 
ized by  a downward  momentum  transfer.  The  Reynolds  stress 
may  be  represented  by 

T = -p  <u'w'>  (8) 

where  u'  = fluctuating  horizontal  wind  velocity 
w’  = fluctuating  vertical  wind  velocity 
p = density  of  air 

It  is  convenient  to  express  Reynolds  stress  in  terms  of  the 
friction  velocity  so  that 

T = pU*^  (9) 

where  Uj.j  is  constant  throughout  a region  of  constant  momen- 
tum flux.  Hence,  is  a measure  of  the  downward  transfer 

of  momentum  in  the  lower  50  meters  of  the  atmosphere.  Over 
the  ocean  an  increase  in  the  near  surface  winds  would  lead 
to  a greater  momentum  and  energy  transfer  for  surface  wave  and 
sea-salt  aerosol  production.  The  relationship  between  the 
*■  jrbulent  transfer  of  heat  and  moisture  in  the  marine  boundary 
. -r  and  the  generation  and  transfer  of  aerosols  is  not  well 
and,  unfortunately,  is  not  investigated  in  this  study. 
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B.  MOMENTUM  TRANSFER,  U,^ , RELATIONS 

A thorough  discussion  of  the  boundary  layer  expressions 
is  presented  in  several  references,  e.g.  Lumley  and  Panofsky 
(1964).  The  similarity  approach  of  Monin  and  Obukhov  (1954) 
is  used  to  define  a representative  length  scale,  L , for  the 
surface  layer  of  the  atmosphere. 


kg 

where  g = gravitational  acceleration 

T = ambient  temperature 
k = von  Karman  constant  = 0.35 

The  selection  of  the  Monin-Obukhov  length  as  a stability  scal- 
ing parameter  is  based  on  the  assumption  that  friction  veloc- 
ity,  U,.j  , and  vertical  heat  flux  (w'T'  ) are  constant  in  the 
surface  layer.  This  scaling  length,  using  dimensional  analy- 
sis, leads  to  the  development  of  a dimensionless  function, 
(})^(Z/L)  , which  can  be  used  to  represent  the  mean  horizontal 
wind  variation  with  height,  du/dZ,  in  the  surface  layer.  The 
following  expression  is  the  empirical  relationship  for  the 
wind  shear  in  this  development, 

dG 

As  vertical  turbulent  heat  flux  (w'!”  ) decreases  to  zero, 
indicating  neutral  stability,  must  approach  1 if 

Equation  (11)  is  to  take  on  its  expected  form  under  neutral 
conditions.  Assuming  that  convective  mixing  is  negligible 
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under  neutral  conditions  it  follows  that  for  values  of  ({) 

m 

(Z/L)  near  1 or  Z <<  1 mechanical  turbulence  is  of  primary 
importance.  Thus,  the  absolute  magnitude  of  L becomes  an 
indicator  of  the  vertical  extent  to  which  mechanical  turbu- 
lence controls  the  turbulent  regime. 

Observational  experiments  by  Businger  et  (1971)  pro- 

duced a definite  relationship  between  the  Richardson  number, 


R. 


goe^/3z) 

5(9u/3Z)^ 


(12) 


and  the  Monin-Obukhov  length,  L , where  0^  is  the  virtual 
temperature.  The  following  expressions  are  approximations 

for  the  unstable  and  stable  conditions  respectively. 


Z/L  = R^  (13) 

R. 

Z/L  = — (14) 

1-aR^ 

where  a is  an  empirically  derived  constant  equal  to  0.5. 

Of  interest  in  this  study  is  the  rate  of  viscous  molecular 
turbulent  kinetic  energy  dissipation,  e . Wyngaard,  e^t  al . 
(1971)  considered  the  dependence  of  e on  momentum  fluxes 
and  height  in  deriving  the  following  empirical  expression 

e = U*^/kZ  (Z/L)  (15) 

Since  Z/L  and  R^  are  functionally  related,  equations  (11) 
and  (15)  can  be  rewritten  as 

du 
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and 


(17) 


e = U*^/kZ  (R^) 

where  and  are  stability  corrections  equal  to  1 

under  neutral  conditions.  In  near  neutral  conditions,  the 
turbulent  kinetic  energy  production  is  assumed  to  be  equal 
to  the  rate  of  molecular  dissipation  of  turbulent  kinetic 
energy  and  from  equations  (16)  and  (17)  the  following  rela- 
tion is  valid 

e = U*^Ou/3Z)  (18)  ■ 

Assuming  neutral  conditions,  the  combinations  of  equations 
(16)  and  (18)  yields 

U*  = (ekZ)^^^  (19) 

Now  the  friction  velocity  U*  can  be  estimated  from  either 
mean  wind  profiles  using  the  integrated  form  of  equation  (16) 
or  from  velocity  fluctuation  data  involving  turbulent  energy 
dissipation  by  using  equation  (19).  The  latter  approach  is 
used  in  this  study. 
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IV.  DATA  COLLECTION 


A.  DURATION  AND  LOCATION 

Aerosol  and  meteorological  data  for  this  study  were 
made  available  through  Calspan  Corporation,  Buffalo,  New 
York,  from  two  separate  experiments.  During  a ten  day  per- 
iod in  February  1977,  Calspan  Corporation  provided  limited 
meteorological  and  cloud  physics  support  during  a study  of 
marine  boundary  layer  phenomena  conducted  on  the  Gulf  of 
Mexico  (Mack  and  Katz,  1977).  The  experiment  was  performed 
on  the  Naval  Coastal  Systems  Laboratory's  (NC3L)  offshore 
platform  "Stage  I"  located  approximately  20  km  SW  of  Pa..ama 
City,  Florida  as  depicted  in  Figure  10. 

A second  experiment  which  provided  d^ta  for  this  study 
was  conducted  along  the  coastal  waters  of  Southern  Califor- 
nia (Figure  11)  during  a 12  day  period  in  July  1977  aboard 
the  Naval  Postgraduate  School  (NPS)  R/V  Acania.  Under  con- 
tract from  NPS,  Calspan  Corporation  provided  limited  meteor- 
ological and  aerosol  physics  support  during  a study  of  air 
quality  parameters  and  marine  boundary  layer  characteristics 
(Mack,  1977).  This  region  contains  primary  shipping  lanes 
and  a number  of  drilling  platforms  all  of  which  contribute 
to  atmospheric  contamination. 

The  following  discussion  will  be  limited  to  equipment 
used  to  measure  the  meteorological  parameters  actually 
analyzed  in  this  study.  A listing  of  the  Panama  City  and 
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Location  of  NCSL  Offshore 


SANTA  BARBARA 
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Figure  11.  Location  of  Southern  California  Cruise 


Southern  California  data  may  be  found  in  Tables  V and  VI , 
respectively,  at  the  end  of  the  text. 


B.  PANAMA  CITY  INSTRUMENTATION 

"Stage  I"  provided  a stable  platform  for  measuring  the 
meteorological  parameters  necessary  to  describe  and  study 
the  aerosol  distribution  and  behavior  in  the  marine  boundary 
layer.  The  instrumentation  installed  by  Calspan  included  a 

1 

Sling  Psychrometer , Bechman-Whitley  wind  system,  Gardner 
small  particle  detector,  and  Royco  Model  225  Particle  Counter. 

The  wind  speed  and  direction  was  monitored  continuously  at 
the  20  meter  level  while  wet  and  dry  bulb  temperatures  were 
obtained  hourly  at  the  17  meter  level.  A Foxboro  temperature 
system  (4  sensors)  provided  continuous  temperature  measure- 
ments at  4 levels;  sea  surface,  4.5,  9.0  and  24.5  meters. 

This  data  was  recorded  in  an  hourly  log.  Ten  minute  averaged 
aerosol  size  spectra  were  obtained  continuously  with  the 
Royce  counter  at  the  17  meter  level,  and  a printout  of 
aerosol  concentration  in  5 size  intervals  was  provided  every 
ten  minutes.  The  Gardner  Counter  measured  the  concentration 
of  particles  greater  than  .0025  y diameter  on  an  hourly 
basis . 

The  majority  of  the  time  the  Royco  instrument  operated 
in  "threshold"  piode  where  number  concentration  (per  2.8  liters) 
of  particles  greater  than  the  following  size  ranges  were 
measured:  0.5  ym,  0.7  ym,  1.4  ym,  3.0  ym,  and  5.0  ym  diameter. 

For  a shorter  period  of  time  the  instrument  was  operated  in 
the  "window"  mode  producing  number  concentrations  between  the 
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above  size  ranges.  The  particle  counter  and  sensor  are  shown 
in  Figure  12.  The  environmental  air  was  drawn  continuously 
through  a sampling  line  of  3 meter  length  and  5 cm  inside 
diameter.  The  flow  rate  through  the  counter's  sensing  volume 
was  set  at  2.8  liters  per  minute. 

The  Royco  Model  225  sampler  utilizes  a near  forward 
scattering  optical  system  (Figure  13)  which  is  ideal  for  moni- 
toring large  volumes  of  ambient  gases  where  suspended  par- 
ticles can  vary  widely  in  composition,  size,  and  optical 
properties.  The  aerosol  is  drawn  through  the  sensor  into  a 
beam  of  focused  light.  As  each  particle  passes  through  the 
illuminated  volume,  it  scatters  a pulse  of  light  which  is 
then  detected  by  a photomultiplier  tube.  The  photomultiplier 
output  is  then  processed  elctronically  to  produce  a pulse 
height  spectrum  from  which  the  particle  size  spectrum  is 
deduced.  The  height  of  each  pulse  is  proportional  to  the 
square  of  the  diameter  of  the  particle. 

Whitby  and  Liu  (1973)  note  that  the  important  character- 
istics of  an  optical  counter  are  the  sampling  flow  rate  and 
the  size  of  the  optical  viewing  volume-.  The  sampling  flow 
rate  determines  the  minimum  counting  period  needed  to  obtain 
a statistically  accurate  count,  and  the  size  of  the  optical 
viewing  volume  determines  the  maximum  aerosol  concentration 
the  instrument  can  accept  without  loss  of  particle  count  due 
to  "coincidence",  i.e.,  the  loss  of  particle  count  due  to 
the  presence  of  more  than  one  particle  in  the  optical  viewing 

volume  at  the  same  time.  The  viewing  volume  of  the  Royco  225 
3 

IS  4.0  mm  and  the  collection  aperture  half  angle  is  25  degrees. 
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Figure  13.  Near  Forv/ard  Scattering  Optical  System 


C-48 


This  model  is  also  equipped  with  a sheath  air  inlet  which 
diverts  part  of  the  aerosol  stream  through  an  external  fil- 
ter before  reentry  to  the  viewing  volume.  This  sheath  im- 
proves the  performance  of  the  instrument  by  preventing  the 
recirculation  of  particles  in  the  optical  chamber  and  by 
confining  the  aerosol  stream  to  a narrower  region.  Thus, 
the  broadening  of  the  pulse  spectrum  due  to  variation  in 
illuminating  intensity  is  reduced. 

C.  SOUTHERN  CALIFORNIA  INSTRUMENTATION 

The  location  of  the  sensors  aboard  the  R/V  Acania  are 
shown  in  Figure  14.  Again,  a Royco  Model  225  Optical  Par- 
ticle Counter  was  used  to  measure  the  aerosol  concentration 
of  the  coastal  marine  boundary  layer.  This  instrument  was 
operated  continuously  in  the  threshold  mode  where  number 
concentration  (per  .28  liters)  of  aerosols  greater  than  the 
following  size  ranges  were  measured:  0.3  pm,  0.6  pm,  1.2  pm, 
3.0  pm,  and  5.0  pm  diameter.  The  mainframe  and  sensor  were 
located  near  the  bridge  of  the  Acania  with  the  origin  of  the 
sampling  line  positioned  forward  of  the  pilot  house  roof  at 
a height  of  7 meters  above  the  sea  surface.  The  sampling 
line  was  6 meters  long  with  an  inside  diameter  of  5 cm.  The 
air  was  sampled  through  the  viewing  volume  at  a rate  of  .28 
liters  per  minute.  A Gardner  small  particle  detector  was 
again  used  to  measure  the  Aitken  nuclei  concentration. 

A sling  psychrometer  was  used  to  measure  the  wet/dry 
bulb  temperatures  and  relative  humidity  determined  from 
psychometric  tables  for  a height  of  5 meters.  The  mean  wind 
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20.5 


measurements  were  obtained  at  four  levels  using  cup  anemo- 
meter wind  profile  register  systems  supplied  by  the  NPS. 
Calspan  recorded  the  wind,  humidity,  and  aerosol  measure- 
ments in  an  hourly  log. 

Velocity  fluctuation  m.easurement s were  obtained  with 
Thermo-Systems  Model  1210  hot  wire  anemometer  probes  mounted 
with  hot  film  sensors  (platinum  coated,  60  mil  quartz  fibers) 
installed  by  the  NPS.  The  anemometer  was  a Thermo-Systems 
Model  1054B.  The  sensors  were  small  enough  to  resolve  the 
viscous  dissipation  scale  without  making  corrections  for 
wire  length.  Wind  fluctuation  data  were  recorded  on  a 14 
channel  tape  recorder.  The  placement  of  these  sensors  re- 
quired exceptionally  long  cable  runs.  Therefore,  adjustments 
were  made  in  the  bridges  for  resistance  and  capacitance 
of  the  wirelength  to  insure  a correct  response. 

The  mean  and  fluctuation  wind  data  were  logged  by  the 
NPS  developed  MIDAS  (Microprogrammable  Integrated  Data 
Acquisition  System) . This  system  is  fully  automated  to 
sample  the  tailored  list  of  sensors  every  30  seconds  and 
20  minute  averaged  output  values  were  printed. 


aO*A054  239  NAVAL  POSTONAOUATC  SCf<K>OL  MONTENEY  CALIF  p/9  13/9 

CXPERIMLNTAL  INVESTX9AT10N  OF  TM£  MARINE  BOUNDARY  LAYER  IN  SUPP— ETC(U) 
FCB  70  9 C SCHACHEA*  C H FAIRALL 


V.  ANALYSES  PROCEDURES 


A.  VELOCITY  FLUCTUATION  ANALYSIS 

The  dissipation  of  turbulent  kinetic  energy,  e , can  be 
related  to  the  mean  wind  velocity  at  any  given  level,  u , 
and  the  RMS  value  of  the  velocity  fluctuation,  u'^  , in  a 
frequency  band  specified  by  a lower  frequency  limit,  fj^  , and 
an  upper  frequency  limit,  f^  (Fairall,  ^ , 1977).  The 

relationship  is 


e 


(4/3)  (u'  ) ^ 

(u/2TT)Cf  f -2/3^3/2 

jC  U 


(20) 


In  this  procedure  recordings  were  made  of  both  the  cup 
anemometer  wind  speed  and  the  corresponding  hot  wire  voltage 
output.  The  sensor  wind  speed  is  given  by 

V z vj  + B(G)’®  (21) 

2 

where  v is  the  hot  wire  voltage  output,  and  and  B 

are  constants  obtained  by  laboratory  calibration  using  a TSI 
Model  1125  Calibrator.  Differentiation  of  equation  (21)  pro- 
duces the  following  relationship  between  the  velocity  fluc- 
tuation and  the  voltage  fluctuation: 


_ 4v(u)^ 


RMS 


RMS 


(22) 
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Substitution  into  equation  (20)  yields 


e 


(u/2TT)Cf  - f -2/3^3/2 

A/  LI 


(23) 


Values  of  = 5 Hz  and  = 200  Hz  were  selected 
for  the  cruise  and  since  amplifiers  with  known  gains,  G , 
were  required,  further  reduction  leads  to 


G = (3.53xlO^)CV^^+  B(u)^3^^^(u)’®[vJy^g/BG]^  (24) 

The  friction  velocity,  U*  , was  then  calculated  from  equa- 
tion (19)  for  each  of  three  levels  and  averaged  to  produce 
over  400  values  from  19-27  July.  Voltage  fluctuation  data 
from  level  3 proved  to  be  erroneous  and  were  not  included 
in  the  calculations.  Obviously  erroneous  values  of  U* 
owing  to  erratic  behavior  were  also  neglected.  values 

were  then  averaged  about  the  aerosol  observation  times  to 
correspond  to  a given  aerosol  distribution. 

B.  AEROSOL  ANALYSIS 

Analyses  were  performed  on  215  aerosol  samples  during 
the  SC  cruise  which  were  confined  to  the  time  period  of  the 
valid  velocity  fluctuation  measurements.  The  observations 
included  date  and  time,  humidity,  relative  wind  speed  and 
direction,  ship's  speed  and  heading,  Aitken  concentration, 
and  aerosol  concentration  as  determined  by  the  Royco  225 
optical  counter  (Table  VI).  Wind  and  ship's  speeds  were  re- 
corded in  knots . 
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The  analyzed  aerosol  observations  for  the  PC  experiment 
were  limited  to  137  cases  during  the  period  18-23  February. 
Cold  frontal  passage  at  approximately  OOOOZ,  24  February  and 
subsequent  advection  of  continental  dust  through  25  February 
were  reasons  for  neglecting  the  aerosol  samples  for  these 
days.  Aerosol  counts  prior  to  18  February  were  determined 
with  the  Roy CO  instrument  in  the  window  mode  and  were  not  in- 
cluded in  this  study.  Observations  were  generally  made 
hourly  and  recorded  in  a log.  They  included  date  and  time, 
humidity,  wind  speed  and  direction  (knots),  Aitken  concentra- 
tion, and  data  from  the  optical  particle  counter  (Table  V). 

Computer  programs  were  developed  to  plot  the  aerosol  size 

distribution  as  a function  of  radius  (R)  in  microns  versus 

- 3 

dN/d  log  R (cm  ) where  N is  the  number  of  particles 
greater  than  a given  radius  as  measured  by  the  Royco  instru- 
ment. The  program  also  included  provisions  to  plot  size  dis- 
tributions predicted  by  Fitzgerald's  model.  For  this  the 
observed  relative  humidity  and  wind  speeds  were  used  with 
equation  (7).  Initially  the  average  aerosol  distributions 
for  both  the  SC  and  PC  experiments  were  computed  and  compared 
to  the  respective  predicted  model  distributions. 

Subsequently,  the  variations  in  the  average  aerosol  dis- 
tributions with  respect  to  four  different  categories  of  wind 
speed,  relative  humidity,  and  friction  velocity  were  plotted 
for  the  SC  data.  The  categories  chosen  for  each  of  the  above 
respective  parameters  are  as  follows:  0-2,  2-5,  5-8,  8-12 
m/sec;  90-99,  80-90,  70-80,  and  60-70  percent;  and  0-.15, 
.15-.  25,  . 25-. 35,  . 35-.  70  m/sec. 
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Friction  velocity  data  was  not  available  from  the  PC 
experiment;  therefore,  variations  in  the  aerosol  distribu- 
tions were  plotted  with  respect  to  categories  of  wind  speed 
and  humidity  only.  Because  of  essentially  different  meteor- 
ological conditions,  the  categories  were  chosen  as  follows: 
0-3,  3-7,  7-10,  and  10-15  m/sec;  and  85-99,  70-85,  55-70  and 
40-55  percent. 

Visual  inspection  of  these  plots  may  indicate  satisfac- 
tory relationships  between  the  aerosol  concentration  and  the 
above  parameters.  However,  a statistical  means  of  viewing 
these  relationships  was  also  deemed  necessary.  Wind  speed, 
humidity,  and  U,-.  values  were  cross  correlated  with  number 
concentration  of  particles  in  graduated  size  ranges.  This 
procedure  was  accomplished  by  a Biomed  Regression/Correla- 
tion computer  program  which  produced  corresponding  correla- 
tion coefficients. 

The  nature  of  the  diurnal  variation  of  the  aerosol  con- 
centration during  the  SC  and  PC  experiments  was  investigated 
in  this  study.  A computer  program  averaged  the  aerosol  con- 
centrations, wind  speeds,  humidities^ and  friction  velocities 
about  each  hour  and  plots  showing  variations  with  time  are 

produced.  The  aerosol  plots  depict  the  number  of  particles 
3 

per  cm  within  specified  size  ranges  versus  time.  The  SC 
data  produced  curves  representing  the  number  of  particles 
between  the  following  size  ranges:  .15-.  30  p,  . 30-.  60  y, 
.60-1.5  y,  and  1.5-2. 5 y radius.  Diurnal  variation  of  con- 
centration for  the  PC  data  utilized  the  following  slightly 
different  size  ranges:  . 25-.  35  y,  . 35-. 70  y,  .70-1.5  y,  and 
1.5-2. 5 y radius. 
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Finally,  diurnal  variations  of  the  aerosol  size  distri- 
bution for  the  SC  and  PC  experiments  were  calculated  using 
techniques  similar  to  those  described  above.  Average  size 
distributions  for  the  following  two  time  periods  were  plotted 
0000-1200  hrs  and  1200-2400  hrs . 

C.  ERROR  ANALYSIS 

The  optical  particle  counter  has  an  advantage  over  the 
membrane  filter  or  impactor  sampling  techniques.  For  example 
the  latter  require  the  samples  to  be  taken  to  a lab  for 
microscopic  inspections  and  the  aerosols  may  possibly  be 
disturbed  or  altered  due  to  contamination.  Although  the 
optical  counter  provides  continuous  "in  situ"  aerosol  measure 
ments,  there  are  ample  causes  for  counting  errors.  Because 
light  scattering  is  a function  of  size,  shape,  and  refractive 
index  of  the  particles,  careful  calibration  is  necessary. 

The  Royco  225  model  counter  used  in  these  experiments 
was  calibrated  using  monodisperse  latex  spheres  of  known  re- 
fractive index  (1.6).  Laboratory  experiments  by  Lieberman 
and  Allen  (1969)  showed  a good  correlation  between  the  theo- 
retical response  curve  for  a near  forward  optical  system  and 
measurements  using  latex  sphere  and  glass  beads  of  refractive 
index  1.6  (Figure  15).  Of  most  significance  is  the  "fold" 
in  the  curve  or  zone  of  multi-valued  response  in  the  region 
of  1 p diameter.  Figure  16  is  provided  to  illustrate  how 
the  response  curve  varies  with  particles  of  different  refrac- 
tive index.  It  is  evident  that  when  measuring  particles  of 
refractive  index  1.6,  a zone  of  ambiguity  exists  between 
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Figure  15»  Theoretical  Response  Curve  and  Experimental 
Results  for  a Forward  Scattering  System 
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RELATIVE  SCATTERING 


Figure  16,  Dependency  of  the  Response  of  a Forward 
Scattering  System  on  Refractive  Index 
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approximately  .3  y and  .6  y radius  and  may  vary  with  the 
aerosol  refractive  index.  Lieberman  and  Allen  (1969)  state 
that  the  instrument  will  still  produce  valid  data  if  the 
zone  is  encompassed  within  a size  range  or  channel.  Since 
the  SC  counter  measured  between  .3  y and  .6  y radius  and 
the  PC  counter  between  .35  y and  .70  y radius,  it  is  assumed 
that  this  multi-valued  zone  is  compensated  for. 

Counting  errors  can  also  arise  from  flow  rate  considera- 
tions. If  the  particle  sizes  are  large  and  the  number  of 
particles  small,  enough  particles  must  be  counted  to  obtain 
good  statistical  resolution.  When  a small  random  number  of 
particles  is  counted,  the  statistical  error  in  counting  is 
equal  to  the  ratio  of  1 over  the  square  root  of  the  number 
of  particles  counted  (Zinky,  1962).  The  counter  should  be 
operated  over  a longer  time  period  (10  minutes)  to  sample 
a larger  volume  or  an  increase  in  the  flow  rate  will  reduce 
the  error.  It  then  seems  quite  possible  that  the  flow  rate 
of  the  counter  used  in  the  SC  cruise  (.28  liter/min)  pro- 
vided too  small  of  a sampling  volume  to  obtain  an  accurate 
count  of  the  larger  particles. 

Zinky  (1962)  also  states  that  a vertically  aligned  inlet 
tube  is  recommended  to  prevent  any  deposition  in  the  line 
due  to  settling.  It  has  already  been  mentioned  that  the 
sampling  lines  used  in  each  experiment  were  considerably 
long  and  aligned  horizontally.  Many  of  the  larger  particles 
may  not  have  remained  airborne  long  enough  to  reach  the 
illuminated  volume. 


Errors  in  the  calculation  of  the  friction  velocity  may 
have  come  from  various  sources.  Since  calculation  of  U* 
from  Equation  (19)  is  only  valid  for  near-neutral  conditions, 
any  substantial  departure  of  the  Richardson  number  from  zero 
would  result  in  inaccurate  values.  The  measurement  of  the 
dissipation  of  turbulent  kinetic  energy  was  large  dependent 
on  the  accuracy  of  the  voltage  output.  The  signal  response 
is  sensitive  to  electromagnetic  energy,  and  any  local  radio 
or  radar  transmission  may  introduce  noise  to  the  system. 
Additionally,  under  the  light  wind  conditions  which  prevailed 
on  the  SC  experiment,  the  lateral  motion  of  the  anemometers 
due  to  ship  pitch  and  roll  may  have  resulted  in  erroneously 
high  readings. 
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VI.  RESULTS 


The  data  from  the  Southern  California  (SC)  cruise  proved 
to  represent  an  atmosphere  somewhat  different  from  a typical 
marine  environment.  The  Aitken  particle  population  averaged 
almost  8500  cm  which  is  about  4 times  higher  than  that  ob- 
served by  Hidy,  et  130  km  west  southwest  of  Los  Angeles. 

This  high  concentration  is  suspected  to  be  due  to  a combina- 
tion of  pollution  from  merchant  ships'  exhaust,  combustion 
from  the  drilling  platforms,  and  offshore  flow  from  the  near- 
by populated  coastal  cities. 

The  average  wind  speed  and  relative  humidity  were  3 m/sec 
and  86  percent,  respectively.  This  data  was  used  to  compute 
the  prediction  from  Fitzgerald's  model  (Eq.  7)  which  is  com- 
pared to  the  average  SC  distribution  in  Figure  17.  The  ver- 
tical bars  represent  one  standard  deviation  either  side  of 
the  mean.  There  is  generally  good  agreement  between  the  two 
below  .4  u radius,  with  a larger  experimentally  observed 
concentration  above  this  range.  Although  sea-salt  production 
should  have  been  minimal  during  this  time  period  because  of 
low  wind  speeds,  the  characteristic  hump  at  around  1 y,  to 
a certain  extent,  reflects  the  contribution  by  sea  salt 
nuclei.  A similarity  exists  here  with  Moore  and  Mason's 
(1954)  observation  of  a discontinuity  where  the  slope  changes 
and  becomes  rather  steep  in  the  region  of  the  larger  size 
range.  The  larger  concentration  in  this  range  may  be  solely 


C-61 


DN  / D(LOG  R)(CM* 


due  to  the  influence  of  atmospheric  contaminants  such  as  com- 
bustion by-products,  soil  dust,  or  smoke.  Considering  pre- 
vious experiments , this  range  does  indeed  contain  a mixture 
of  both  continental  and  marine  aerosols  possibly  resulting 
in  the  increase  over  Fitzgerald's  model. 

As  previously  mentioned,  the  low  flow  rate  of  the  optical 
counter  may  account  for  the  low  concentrations  at  2 y . How- 
ever, since  the  wind  speed  reached  8 m sec  ^ only  6 times, 
this  may  have  been  a truly  representative  concentration  of 
droplets  as  agreement  is  also  shown  with  Fitzgerald's  curve. 

Figure  18  presents  the  synoptic  situation  for  three  days 
at  the  beginning,  middle,  and  end  of  the  experiment.  A per- 
sistent thermal  low  is  located  in  the  desert  area  of 
Southern  California  and  the  isobaric  pattern  off  the  coast 
reflects  a rather  weak  gradient.  Therefore,  smaller  scale 
circulations  should  prevail  in  this  area  of  little  or  no 
synoptic  forcing. 

Plots  showing  the  variations  of  the  average  distribu- 
tions with  respect  to  wind  speed,  relative  humidity,  and 
friction  velocity  (U.;.),  are  shown  separately  in  Figures  19, 
20,  21.  The  number  of  observations  in  each  category  is 
placed  in  parentheses.  These  figures  indicate  that  the  size 
distribution  has  a better  relationship  with  the  relative 
humidity  than  to  the  wind  speed  and  U*  . Correlation  co- 
efficients between  these  parameters  and  the  number  concen- 
tration of  particles  in  a given  size  interval  are  produced 
in  Table  II.  Since  diurnal  variations  tend  to  reflect  a 
negative  relationship  between  relative  humidity  and  wind 
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19  JUL 


23  JUL 


26  JUL 

Figure  18,  Synoptic  Situation  during  the  SC  Experiment 
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Figure  19,  Variation  of  SC  Size  Distribution  with 
Wind  Speed 
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Figure  20,  Variation  of  SC  Size  Distribution  with 
Relative  Humidity 


speed,  the  results  in  this  table  are  not  surprising.  The 
negative  correlation  of  relative  humidity  with  the  concen- 
trations in  the  large  size  range  indicates  that  sedimenta- 
tion of  large  droplets,  which  grow  with  increasing  humidity, 
is  most  important  when  the  wind  speed  and  sea  surface  pro- 
duction of  salt  nuclei  are  weak.  Although  these  larger  drop- 
lets also  exhibit  a small  positive  correlation  with  U,.5 
while  the  wind  speed  correlation  remains  negative,  this  re- 
sult does  not  appear  to  be  significant. 

An  attempt  was  made  to  examine  the  influence  of  stability 
on  the  size  distribution.  The  summer  months  are  character- 
ized by  the  occurrence  of  stratus  and  fog  off-shore  below 
the  marine  inversion.  Two  days  are  compared  with  the  assump- 
tion that  they  represent  the  unstable  and  stable  atmospheres. 
According  to  the  daily  observation  log,  stratus  clouds  in 
the  morning  becoming  stratocumulus  by  afternoon  were  observed 
on  19  July.  26  July  was  characterized  by  clear  skies.  The 
average  distribution  for  both  days  is  presented  in  Figure  22. 
The  correlation  coefficients  between  concentration  and  wind 
speed  and  U*  show  a trend  toward  positive  values  from  the 
stable  to  the  unstable  day  with  U*  eventually  becoming 
positively  correlated  in  the  unstable  day  (Table  III).  The 
increase  in  the  size  distribution  on  26  July  in  the  size 
range  greater  than  .3  y seems  to  be  due  to  increase  in  the 
average  wind  speed  and  occasional  gustiness  as  whitecaps  were 
reported  during  the  afternoon.  The  stable  stratification 
assumed  in  this  case  allows  generated  sea-salt  nuclei  to 
accumulate  and  the  concentration  to  increase  at  the  7 meter 
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Figure  22.  Average  Size  Distributions  on  19  July  and 
26  July 
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Table  HI.  Correlation  Coefficients  for  19  July  and 
26  July 
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level.  The  lower  average  wind  speed  associated  with  the  un- 
stable period  does  not  allow  for  much  sea-salt  production. 

An  unstable  atmosphere  can  lead  to  convective  processes 
which  may  vertically  transport  aerosols  and  create  higher 
concentrations  at  an  upper  level  as  proposed  by  Ericksson 
(1959)  and  Toba  (1965a  6 b).  Hence,  a decrease  in  the  size 
distribution  on  19  July  is  observed.  This  evidence  gives 
credence  to  the  possibility  that  friction  velocity  is  a 
better  indication  of  aerosol  size  distribution  than  wind 
speed.  On  both  days  the  correlation  of  the  concentration 
with  humidity  is  lowest  in  the  largest  size  interval.  This 
relationship  is  most  pronounced  on  the  unstable  days  and  may 
be  explained  by  sedimentation  due  to  mixing  and  resulting 
increased  coalescence. 

The  averaged  diurnal  variations  of  wind  speed,  relative 
humidity,  friction  velocity,  and  aerosol  concentration  are 
shown  in  Figures  23,  24,  and  25.  Again  the  negative  rela- 
tionship of  aerosols  to  wind  speed  and  U*  in  the  size 
range  of  generally  less  than  1 y is  indicated.  A satisfac- 
tory relationship  with  relative  humidity  is  not  evident  and 
this  is  probably  due  to  transport  by  a secondary  circula- 
tion. A land-sea  breeze  type  of  effect  could  account  for 
the  observed  decrease  .in  concentration  of  the  particles 
"^smaller  than  1.5  y.  As  the  heating  over  the  land  generates 
an  on-shore  flow  along  the  coast,  the  wind  increases  and 
persists  through  the  afternoon.  The  average  wind  direction 
derived  from  the  observations  of  five  random  days  during 
the  experiment  is  shown  in  Figure  26.  A westerly  wind  is 


C-72 


FRICTION  VELOCITY  (m/s) 


6 13  18  34 

TIME 


Figure  24.  Average  SC  Diurnal  Variation  of  Friction 
Velocity 
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seen  to  dominate  during  the  peak  wind  periods.  The  decrease 
in  the  aerosol  population  may  be  explained  by  a horizontal 
divergence  effect  in  the  marine  boundary  layer  as  largest 
accelerations  are  found  near  the  coast.  The  average  size 
distributions  displayed  in  Figure  27  reflect  the  decrease  in 
the  aerosol  population  due  to  this  sub-synoptic  circulation. 
Although  the  relative  humidity  increases  slightly  in  the 
early  evening  hours,  the  smaller  nuclei  show  a stronger 
relationship  with  the  wind  speed.  This  again  implies  that 
a large  part  of  the  coastal  marine  aerosol  is  of  continental 
origin.  The  outflow  of  circulation  aloft  is  probably  respon- 
sible for  the  introduction  of  continental  particulates  to 
the  marine  environment.  The  minor  peaks  in  the  small  particle 
concentration  and  also  the  somewhat  greater  increase  in  the 
large  particles  during  the  afternoon  should  be  attributed 
to  sea-salt  production. 

The  Panama  City  (PC)  observations  more  closely  resembled 
a marine  environment.  The  Aitken  particle  count  was  lower 
and  averaged  2600  cm  while  the  distribution  curve  showed 
a marked  change  from  the  Southern  California  data.  Winter 
time  synoptic  scale  features  predominate  in  this  region  of 
the  Gulf  Coast.  Cold  frontal  passages  and  an  accompanying 
influx  of  continental  air  into  the  Gulf  of  Mexico  are  fre- 
quent occurrences.  Subsequent  movement  of  the  high  pressure 
ridge  into  Florida  and  of f its  eastern  seaboard  provides  the 
circulation  which  reestablishes  moist  southerly  flow  and 
return  of  the  marine  aerosol.  Figure  28  provides  the  synop- 
tic analyses  for  the  period  of  the  experiment.  Stable 
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Figure  28,  Synoptic  Situation  during  the  PC  Experiment 
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conditions  prevailed  in  the  Panama  City  area  at  the  begin- 
ning of  the  period;  but,  after  frontal  passage  early  on 
20  February,  south  to  southeasterly  flow  developed  and  per- 
sisted for  the  remainder  of  the  experiment.  The  influx  of 
this  warm,  moist  air  contributed  a destabilizing  effect  in 
the  lower  levels  of  the  marine  boundary  layer. 

The  average  wind  speed  and  relative  humidity  for  PC 
were  8.4  m/sec  and  71  percent,  respectively.  Fitzgerald’s 
curve  for  these  average  conditions  and  the  average  aerosol 
distribution  for  the  entire  period  are  shown  in  Figure  29. 
Good  agreement  exists  only  for  particle  size  range  greater 
than  .9  p radius.  The  observed  concentrations  are  approxi- 
mately an  order  of  magnitude  lower  than  Fitzgerald's  pre- 
diction for  aerosols  smaller  than  .5  p radius.  A signifi- 
cant aspect  of  the  distribution  is  the  positive  slope 
observed  between  approximately  .5  p and  1 p radius  which 
appears  to  be  the  result  of  sea-salt  production.  Actually, 
good  agreement  is  shown  with  Blifford's  (1970)  observation 
off  the  Pacific  coast  with  respect  to  both  slope  and  number 
concentration . 

Plots  showing  the  effect  on  the  average  distributions 
due  to  wind  speed  and  relative  humidity  separately  are  shown 
in  Figures  30  and  31.  Relatively  good  correlations  seem  to 
exist  between  these  parameters  and  the  aerosol  distributions. 
Correlation  coefficients  are  presented  in  Table  IV.  The 
synoptic  scale  effects  predominate  over  diurnal  variations 
and  wind  speed  and  relative  humidity  are  both  positively 
correlated  to  the  concentration.  The  highest  correlation  of 
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Figure  31,  Variation  of  PC  Size  Distribution  with 
a. Relative  Humidity 
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concentration  to  humidity  at  Panama  City  is  witnessed  in  the 
1.5  p to  2.5  u interval.  This  result  may  indicate  that 
equilibrium  tends  to  exist  between  production  and  sedimenta- 
tion in  this  interval  as  hypothesized  by  Iloore  and  Mason 
(1954-).  Disagreement  exists  in  that  the  steeper  negative 
slope  is  found  during  periods  of  strongest  wind.  The  plot 
showing  the  effect  of  relative  humidity  on  the  size  distri- 
butions results  in  small  variations  in  the  .25  p - .35  p 
interval.  This  probably  indicates  that  the  majority  of  these 
particles  represent  a mixture  of  continental  and  marine 
nuclei . 

The  stability  influence  was  investigated  by  comparing 
observations  on  18  February  and  21  February.  Temperature 
measurements  at  various  levels  on  the  platform  made  it  pos- 
sible to  examine  the  lapse  rates  and  determine  the  stability. 
The  average  distributions  for  these  days  grouped  according 
to  wind  speed,  and  respective  correlation  coefficients  are 
shown  in  Figures  32  and  33.  The  low  humidities  on  21 
February  resulted  from  the  earlier  intrusion  of  continental 
air,  but  southeast  to  southwest  flow  persisted  most  of  the 
day.  Although  this  trajectory  helped  to  advect  in  warmer 
air,  production  of  sea-salt  dropped  off  as  the  wind  decreased 
considerably  below  7 m/sec.  A much  larger  decrease  is 
observed  in  the  distribution  curve  on  21  February  as  compared 
to  18  February  when  the  wind  speed  decreased  below  7 m/sec. 
This  agrees  well  with  Moore's  (1952)  finding  that  the  change 
in  opacity  is  well  marked  during  periods  of  low  humidity. 

Also  the  decrease  in  the  slope  of  the  curve  between  .5  p and 
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Figure  32,  Correlation  Coefficients  and  Variation  of  the 
Size  Dlstrihution  with  Wind  Speed,  18  February 
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Figure  33,  Correlation  Coefficients  and  Variation  of  the 
Size  Distribution  with  Wind  Speed,  21  February 
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1 p appears  to  be  a function  of  decreased  relative  humidity. 
This  again  reinforces  the  premise  that  sea-salt  nuclei  pre- 
dominate in  the  size  range  above  approximately  .7  p. 

The  correlation  between  concentration  and  wind  speed 
for  all  size  ranges  is  greatest  on  the  unstable  day.  As 
with  the  SC  data,  this  is  consistent  with  momentum  and  diffu- 
sion theory.  Again  the  small  correlation  with  relative 
humidity  exhibited  by  the  larger  nuclei  is  probable  caused 
by  growth  and  sedimentation  in  the  absence  of  significant 
sea-salt  production.  When  generation  was  occurring  on  18 
February,  the  large  particles  exhibited  the  largest  corre- 
lations with  humidity  and  wind  speed.  This  stable  stratifi- 
cation evidently  was  produced  by  a previous  frontal  passage 
and  northerly  flow  of  cold  air  and  accompanying  continental 
particulates.  Therefore,  a large  portion  of  the  aerosol  at 
the  beginning  of  the  experiment  may  have  been  composed  of 
non-hygroscopic  material. 

Figures  34  and  35  display  the  average  diurnal  changes 
in  wind  speed,  relative  humidity,  and  aerosol  concentration. 
Again  positive  correlations  are  noted  as  relative  humidity 
and  wind  speed,  although  containing  quite  a bit  of  scatter, 
tend  to  vary  accordingly.  Of  most  significance  would  be  the 
obviously  high  concentration  of  droplets  in  the  .7  p - 1.5  p 
range.  Noting  that  the  average  wind  seldom  went  below  7 m/ 
sec,  this  would  indicate  that  sea-salt  nuclei  production 
is  greatest  in  this  size  range.  A diurnal  representation 
of  the  average  aerosol  distribution  is  presented  in  Figure  36. 
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Figure  35,  Average  PC  Diurnal  Variations  of  Particle 
Concentrations 
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Figure  36,  Average  PC  Diuj.nal  Variation  of  the  Aerosol 
Size  Distribution 
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Any  transport  of  aerosols  due  to  a land-sea  breeze  effect 
should  be  ruled  out  as  a satisfactory  relationship  does  not 
seem  to  exist. 
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VII.  CONCLUSIONS 


The  coastal  marine  aerosol  is  shown  to  be  a highly  vari- 
able function  of  the  interaction  between  synoptic  and  meso- 
scale  processes.  Important  meteorological  parameters  such 
as  wind  speed,  relative  humidity,  and  stability  are  dependent 
upon  secondary  circulations  between  land  and  sea  in  the 
absence  of  large  scale  forcing. 

The  minimum  concentration  in  the  size  distribution  curves 
at  .4  y - .5  p radius  may  indicate  that  this  size  is  indeed 
the  transition  zone  between  the  two  bubble  bursting  sea-salt 
producing  mechanisms.  Since  the  slope  on  either  side  of 
this  zone  is  steeper  during  the  Panama  City  experiment,  wind 
speeds  of  greater  than  7 m/sec  result  in  the  generation  of 
sea-salt  particles  larger  than  .25  y radius.  Sedimentation 
of  particles  larger  than  1.5  y appears  to  be  most  significant 
during  periods  of  low  wind  speed.  During  strong  winds  a 
state  of  equilibrium  between  sedimentation  and  production 
exists  for  these  larger  particles. 

Relative  humidity  variations  have  the  largest  effect  on 
the  aerosol  size  distribution  in  the  absence  of  sea-salt 
production.  The  concentration  of  the  coastal  marine  aerosol 
is  most  sensitive  to  wind  speed  effects  at  low  relative  hu- 
midity. Friction  velocity  seems  to  be  a better  indication 
of  the  aerosol  size  distribution  than  wind  speed  under  un- 
stable atmospheric  conditions.  Also  during  light  wind 
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periods,  instability  appears  to  result  in  a decrease  in 
concentration  at  the  observation  height.  Enhanced  diffu- 
sion during  periods  of  sea-salt  production  causes  vertical 
transport  of  sea-salt  from  the  sea  surface  and  an  increase 
in  concentration. 

Any  effect  of  surface  organic  film  possibly  suppressing 
the  production  of  small  sea-salt  particles  could  not  be 
examined  because  of  the  absence  of  significant  generation 
off  the  Southern  California  coast. 
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Table  VI.  Southern  California  Data 
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